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1. INTRODUCTION 



In a recent article Foutz (1980) introduced a new test for goodness-of- 
fit, to be called the test in the sequel. Although the test was proposed 
for fitting a continuous p-variate distribution, it applies equally well to 
univariate problems. The null distribution of the test statistic was shown 
to be distribution free as well as being independent of the number of vari- 
ates p. Foutz obtained an integral representation for the null CDF of F ; 

n 

explicit expressions for this CDF were given for sample sizes 2 and 3. 

Closed form solutions for the CDF for larger sample sizes are quite hard 

to derive and Foutz has provided a large sample normal approximation to the 

null distribution of F . 

n 

In a preliminary comparison with ten replications of 50 simulated samples 
from (I) a mixture of uniform distributions and (Li) a standard normal distri- 
bution, Foutz found that the F test outperformed both the Chi squared test 

n 

and the Kolmogorov-Smirnov (K-S) test. 

In this paper we present the results of an extensive investigation to 

compare the three goodness-of-f it tests, the Chi squared test, the K-S test 

and the F test. Members from three families of distributions, viz., the 
n 

family of asymmetric stable dis tributions, mixtures of normals, and the 
Pearson family have been selected to represent the true underlying distri- 
bution of the samples. T.he goodness-of-f it tests are applied to test the 
hypothesis that the samples are from a standard normal distribution. The 
measure of comparison used in the study is the empirical power, based on 
5000 replications of each of the tests. 

The Chi squared and Kolmogorov-Smirnov statistics were computed using 
the IMSL + routines GFIT and NKS1. The number of cells used in the Chi 

"^International Mathematical and Statistical Libraries, 7500 Bellaire 
Boulevard, Houston, TX 77036 
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squared computation was 6, 8, and 8 for sample sizes of 20, 30, and 50, 
respectively. 

A brief discussion of the F^ test is given in Section 2 and a descrip- 
tion of the simulation is in Section 3. The results of the simulation are 
presented in Section 4. FORTRAN codes used for the simulation and detailed 
tables of simulation results are in Appendices I and II. 
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2. F TEST 
n 

The F^ test is based on a comparison of a continuous empirical distri- 
bution function (CEDF) with the hypothesized CDF. The CEDF is obtained by 
,f spreading M the total mass over "statistically equivalent blocks" generated 
by the sample. As shown in Anderson (1966) and Foutz (1980), given the 
order statistics of a random sample of size (n-1) , n statistically equiva- 
lent blocks that partition the sample space can be constructed in many 
different ways by choosing what are called cutting functions. An intuitively 
appealing set of blocks, which is the one used in this study is obtained by 
choosing the identity functions for the cutting functions. In this case, 
the n statistically equivalent blocks are = (-<», x^], ^ = (x^^x^jl,.. 

B = (x, n ,°°) where x /ix , , . . . , x. are the order statistics of 

n (n-1; (1; (2; (n-1; 

a sample of size n-1. The CEDF is constructed by spreading a mass — cont- 

n 

tinuously and in the same proportion as the hypothesized CDF over each 

block. If H is the hypothesized CDF and H the CEDF, the test statistic 
o n 

F is defined as 
n 

F = Sup |h fx) - H (x) | (1) 

n x r 1 n o 1 

Let D^, i = 1,2, ... ,n be the probability contents of the blocks B i 

under the null CDF H , i.e., D. = P[x e B Ih ]. A computationally convenient 

o 1 l o 

form for F can be shown to be 
n 

n 1 

F = E max(0, D.) . (2) 

n . . n l 

i=l 

Foutz has provided the null distribution of F in integral form and derived 

n 

closed form solutions for n = 3,4. Simplifying the expressions given by 
Foutz for n = 3,4 yields 



6x 

P(F 3 < x) = j 1 - 3(j - X ) 2 



0<x< — 

J < x < 2/3 
x > 2/3 



(3) 



and 



20x J 



P(F 4 < x ) 



3 2 9 1 

-20x + 18x - fx + tt 
4 16 

1 - 4(f - x) 3 



0 < x < — 

- -4 

1 1 
4 - X - 2 

*>!• 



(4) 



We have also obtained the expression for the null CDF for n = 5, 

4 



P(F 5 < x) =< 



70x 

-105x 4 + 80x 2 - 12x 2 + -j|-x - Y 25 - 
45x 4 . 80x 3 + 2|8 x 2_176 + _|I 

1 - 5(j - x) 4 



0 < x < — 

1 2 

5 < X * 5 

2 3 

5 X - 5 

3 4 

5 < X - 5 

4 

X > 5 * 



(5) 



As is evident the exact distribution is quite difficult to obtain for 
higher sample sizes. A large sample normal approximation, due to Foutz is 
given by 



lim P 
n-*” 



F < x 

n - 



= $ 



n(x - e 



(2e 



-2 % 
- 5e Z ) 2 



( 6 ) 



where $ is the standard normal CDF. For n-1 = 20,30,50 we used this approxi- 
mation to test the hypothesis that a simulated sample from U[0,1], a uniform 
distribution on [0,1], is in fact from that distribution. In 80,000 repli- 
cations, the observed significance level (number of hypothesis rejections/ 
80,000) was consistently smaller than the nominal value as can be seen from 



8 



TABLE 1 



EMPIRICAL SIGNIFICANCE LEVEL OF 
FOUTZ F n TEST USING ASYMPTOTIC APPROXIMATION 
(80,000 REPLICATIONS) 



Sample Size 


20 


30 


50 


Significance Level 








.10 


.0757 


.0800 


.0859 


.05 


.0372 


.0399 


.0428 


.01 


.0082 


.0083 


.0093 
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Table 1. We therefore constructed a Monte Carlo CDF of F for n-1 = 2,3,4, 

n 9 

20,30,50 based on 25,000 computer generated F values; these represent values 

n 

of the F^ statistic for testing the hypothesis that a set of samples from 

U[0,1] is in fact from U [ 0 , 1 ] . A comparison of the Monte Carlo CDF with 

the exact CDF for n-1 = 2,3,4 is provided in Table 2. It can be seen that 

the Monte Carlo CDF provides a reasonable approximation even for small n. 

The power properties of the F^ test detailed in this paper are based on 

the Monte Carlo CDF of F . Critical values obtained from the Monte Carlo 

n 

simulation for significance levels .01, .05, .1 and n-1 = 20,30,50 are in 
Table 3. In Table 4 we present the observed significance level in 225,000 
replications when testing if a set of samples from U[0,1] is in fact from 
U[0,1]. 
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TABLE 2 



MONTE CARLO SIMULATION VS EXACT VALUES OF CDF 





n=3 




n=4 




n=5 






MC 


EXACT 


MC 


EXACT 


MC 


EXACT 


X 














.40 


.7883 


.7867 


.7624 


.7625 


.7607 


.7600 


.45 


.8604 


.8592 


.8725 


.8725 


.8671 


.8693 


.50 


.9182 


.9167 


.9396 


.9375 


.9405 


.9405 


.55 


.9598 


.9592 


.9694 


.9680 


.9768 


.9778 


.60 


.9864 


.9867 


.9860 


.9865 


.9915 


.9920 


.65 


.9989 


.9991 


.9952 


.9960 


.9975 


.9975 
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TABLE 3 



CRITICAL VALUES FOR F n TEST 
OBTAINED BY MONTE CARLO SIMULATION 



Sample Size 


20 


30 


50 


Significance Level 








.10 


.42714 


.41903 


.40816 


.05 


.44865 


.43553 


.42116 


.01 


.48659 


.46579 


.44487 
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TABLE 4 



EMPIRICAL SIGNIFICANCE LEVEL OF 
FOUTZ F n TEST USING MONTE CARLO APPROXIMATION 
(225,000 REPLICATIONS) 



Sample Size 


20 


30 


50 


Significance Level 








.10 


.1006 


.0970 


.1003 


.05 


.0486 


.0486 


.0498 


.01 


.0103 


.0101 


.0102 



13 



3. DESCRIPTION OF SIMULATION 



In our simulation we generated deviates from three families of distri- 
butions. The family of asymmetric stable distributions has been previously 
used by Saniga and Miles (1979) to investigate the power of several goodness- 
of-fit tests. We used the same set of parameter values, a = 1.0(. 3)1.9 and 
B = 0(. 25)1.00 they used. Mixed normal distributions have often been used 
for such tests, and we have included a family which is a composite of N(0,1) 
and N(0,a), a ^ 1, and another set which is a composite of N(.5,l) and 
N(0,a). Pearson distributions considered include a variety of shapes, from 
'hear” normal to U and J shaped. Discussion of the procedures used to 
generate pseudorandom deviates from each of these families follows. 

The sample data was obtained by starting with one or more uniformly dis- 
tributed pseudorandom deviates. These were generated using the IMSL sub- 
routine GGUBS, the basis of which is discussed in Lewis, Goodman, and 
Miller (1969) . 

3.1 Generation of Asymmetric Stable Deviates (Random Stabilized Standard 
Form) 

This family of distributions contains as a special case the normal 
distribution (a = 2 , B = 0) and the Cauchy distribution (a = 1, B = 0) . For 
a < 2 the distributions have infinite variance, which makes them useful for 
determining the ability of a goodness-of-f it test to detect heavy tailed 
distributions. The deviates were generated using the program given in 
Chambers, Mallows, and Stuck (1976). This subroutine, RSTAB, used one deviate 
from U [0 , 1 ] and one exponentially distributed deviate to generate one RSSF 
deviate. 
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3.2 Generation of Mixed Normal Deviates 



Mixed normals of the form (l-y)N(y^ ,a^) + y NCy^jO^) were generated using 

the IMSL routine MDNRIS to convert uniform samples into standard normal 

samples. To obtain a set of N mixed normal variates we proceeded as follows: 

(i) 2N uniform random variates {u^} were generated using GGUBSr (ii) for 

each i = 1,...,N, MDNRIS was used to convert u. to a standard normal z.. 

1 1 

was then transformed to a normal with mean y^ and variance a ^ , or with 

mean y n and variance a n , depending on whether u. fx > Y ■> or u. 1XT < Y, 

2 2 x+N - x+N 

respectively. 



3.3 Generation of Pearson Type I, II Deviates 

The generation of samples from Pearson Type I and II distributions was 
done via table look up and linear interpolation on the inverse cumulative 
distribution function. Sufficient entries to assure four significant decimal 
places in the final answer was achieved adaptively using numerical integration. 
Before discussion of the precise details of the process, we digress for a 
discussion of the Pearson distributions, and particularly types I and II. 

Following Johnson and Kotz (1970), the Pearson probability density function 
p(x) is given by 



1 d£ 
p dx 



a + x 



2 

c + c,x + c 0 x 
o 1 2 



It can be shown that a, c , X- , and c n can be expressed in terms of non- 

o 1 l 

negative parameters 8 , and t ^ ie variance y^, obtaining 

C o = (4B 2 ~ 36 1 )(10e 2 - 126 1 - 18) -1 u 2 

a = c x =/f^ (8 2 + 3)(108 2 - 12^ - 18)' 1 /uj” (7) 

c 2 = (28 2 -36 -l -6)(108 2 - 12^ - 18) _1 . 
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Type I distributions are characterized by 



12 -I 

K: = 4 c i ( C 0 C 2> < 0 » 

while Type II have k = 0 with 6^ = 0, < 3. For Type I and II the prob- 

ability density function is of the form 



m l m 2 

p(x) = K(x - a 1 > X (a 2 - x) Z 



(8) 



Generation of a sequence of pseudorandom Pearson Type I and II deviates 
involves the following steps: (i) generate a sequence of deviates {u_^}, 

from U[0,1]; (ii) transform these to Pearson deviates by finding v. so that 

t T i 

/ p(x)dx = u_^. This necessitates being able to efficiently obtain the 



■/ 



p(x)dx, then we need v_^ = F ^(u_^). We 



inverse CDF, i.e., if F(v) = 

will denote F ^ by G in order to simplify notation. 

Representation of the inverse CDF, G(s), was achieved by linear inter- 
polation in a table generated by numerical integration of p(x). The adap- 
tive process used to assure four decimal place accuracy, i.e., magnitude of 

-4 

the error less than .5 x 10 is described now. 

The error in linear interpolation between points (s^,G_^) and (s^_^, ^i+l^ 

12 M 
is no more than —(s. + ^ - s^) Q 9 , where G 2 = max |g (s) | . 

s . <s<s _ 

1 - - 1+1 

Since the intervals will be small, we can approximate 

and then using the fact that 



G 2 by 



G '(s i+ 1> - G -(s.) 



S..T - S. 
1+1 1 



G'(s) = (F 1 (s))’ = rrr 1 



F'(G) P(G) 



, we obtain 
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p«W p<V 

s i+l - s i 



Thus, an error estimate in (s^, is given by 

1 1 



p ^ G i+l^ P( G p 



s . - s . 

1+1 1 



(S H-1 - s i>' 

8 



p( G i+ l) p( G i^ 



(s 



i+1 



S i } 



(9) 



Potential problems occur if p(G) = 0 , as may happen at and a^. Since we 

— 

require error less than .5 x 10 , it is clear this must be the case when 

i i —4 

- Gj<.5 x 10 ; hence if p(G) is very small in the interval 

we have an alternative scheme for accepting an interval, one which came into 
play for ranges where p is small. 

These ideas were the basis for adaptive construction of a suitable table 



(s^,G^) for the inverse cumulative distribution function. We have 

s = 0, G = a. . We describe the general scheme for obtaining (s . , . , G.,-) 
o o 1 i+l i+l 

given and note how an estimate of G_^ + ^ is generated afterwards. Given 

an estimate for G.,-, the value of As. = s. (1 - s. is obtained by numerical 
i+l l i+l l 



numerical integration of 



/ 



G i + i 



p(x)dx. This is accomplished by the 



adaptive quadrature routine, DCADRE, from the IMSL library. If m^ or m^ are 

negative, subtracting out the singularity was used for intervals near a^ and 

a^, respectively. The double precision version of DCADRE is used and an 

— 6 

absolute error tolerance of 10 is requested. The routine returns As^ and 

-4 

the error estimate (9) is computed. If it is less than .5 x 10 , the 



result is accepted, we set s . , „ = s. + As . , and proceed to the next interval. 

i+l l i 

-4 2 

Suppose the error estimate, E > . 5 x 10 . Since E = 0((AG) ), and 

est est 
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O(AG^) = O(As^), we obtain a new estimate for AG_^ by taking it to be 




.25 x 10 



est 



-4 



+ 1 



This takes the error based on the new AG_^ to approximately midway between its 

— /j. 

current value and .25 x 10 . More than one correction of this sort may be 

required, in particular initially where a reasonable estimate of AG is not 

o 

available. 

Once an interval has been accepted (or rather, a point (s^Gj) we 
increment the interval counter i, and then estimate the new interval size 



AG. by AG i-l 
1 2 



.475 x 10 



est 



-4 



+ 1 1. This yields, based on the above assumptions, 



a AG^ which should give an error for the next interval which is midway between 

—4 

the previous one and .475 x 10 

An initial value for AG is required. Although not very sophisticated, 

o 

a 2 ” a i 

we simply take A G q = — > and depend on the adaptive machinery described 

above to decrease it to meet the error tolerance, or increase successive 
intervals as required for efficient representation. 

It is true that the final value of s^, call it s^, should be equal to one. 
Because of the numerical integration, this is never exactly the case. In 
the worst case, 3^ = .01, = 1-9, we obtained the final value 



s._ - .9999999845 an error of about 1.55 x 10 
N 



-8 



In order to avoid any prob- 
lems due to the table not covering [0,1] exactly the simplest procedure 



was to replace each computed s_^ by s ^/s , thus distributing the error over 

the entire interval and yielding a consistent table. Note that this is 

-4 

well within the error tolerance of .5 x 10 

The procedure has been thoroughly tested for its efficiency in repre- 
senting the inverse CDF as well as for accuracy. For the particular cases 
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in which we were interested, the inverse CDF was represented by a table with 
no more than 729 entries. This occurred for the case 6^ = .25, $2 ~ 3.2, 
which has an inverse CDF with very large slopes. For more gently sloping 
inverse CDF f s we were able to use as few as 101 intervals, as in the case 

= .01, £$2 “ 1*75. Most intervals had an error estimate of between .4 x 10 

-4 

and .5 x 10 , which shows that the interval sizing process we have used 

worked quite efficiently, with few initial estimates being rejected for being 
too large, without also resulting in intervals much too small. 

The routine was checked against the published tables of Johnson, Nixon, 
Amos and Pearson (1963) for many values of the parameters and at most of the 
percentage points. With a few exceptions, where a difference of one in the 
fourth decimal place was noted, the results check exactly. Generally in 
these cases the fifth place was four or five so that the actual error was 
probably well within our tolerance. 



\ 
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4 . RESULTS 



The results of the simulation are summarized in Tables 5-12. The empi- 
rical power, in 5000 replications, of the Foutz F test compared with that 

n 

for the Chi squared test or the K-S test is presented as a percent improve- 
ment in the probability of rejecting the null hypothesis that the distri- 
bution of the samples is the standard normal. A negative entry means that 
the power of the F^ test was smaller than that for the Chi squared test or 
the K-S test, whichever is appropriate. 

The simulation has revealed that the F test is better than the Chi 

n 

squared test which in turn is better than the K-S test when the true distri- 
bution of the samples is heavy tailed. Many such distributions are included 
in the asymmetric stable family as well as the family of mixtures of normals. 
For the mixed normal family, if the two normals involved in the mixture 
differ in the means the K-S test performed better than the F^ test even when 
the variances differed. We now discuss in more detail the results for each 
of the three families of distributions. 

4.1 Asymmetric Stable Family 

The results for the F test versus the Chi squared test are summarized 

n 

in Table 5. The F test outperformed the Chi squared test for n equal to 
n 

21 and 31. When n ■ 51, as a + B increased the performance of the F^ test 
deteriorated as can be seen from the lower right part of Table 5. Another 
general observation is that as the significance level is decreased the improve- 
ment in power for the F^ test is accentuated. 

The comparative figures for the F^ test versus the K-S test are pre- 
sented in Table 6. Here again the F test did much better than the K-S test; 
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the results also indicate that for the asymmetric stable family the Chi 
squared test has a higher power than the K-S test. 

4.2 Mixtures of Normals 

The mixed normal distributions that we considered were of two basic 

types. The first type is of the form (l-y)N(0,l) + yN(0,cr) with a = 2,3,4 

and y = .1, .2, .3, 1.0; note that when y = 1.0 the distribution is not a 

true mixture but N(0,a). The second type is a mixture of the form 

(l-y)N(.5,l) + yN(0,a) with a = 3, y = .2, .3 and a = 4, y = .2. 

The test did significantly better than the Chi squared test except 

for n = 51 and y = 1.0 (see Table 7); in the latter case the Chi squared 

test turned out to be the better of the two tests. 

When the F test was compared with the K-S test the F test turned out 
n n 

to be consistently better (Table 8) . 

In the case of a mixture of the second type, which included a location 

shift (Tables 9, 10), the K-S test proved to be superior to the F test 

n 

while the comparison between the F test and the Chi squared test appeared 

n 

to be inconclusive. 

4.3 Pearson Family 

We chose ten distributions of types I and II to encompass a variety 
of shapes as shown in figure 1; the standard normal is superimposed as a 
dotted curve in each of the graphs. The comparison of the F^ test versus 
the Chi squared test proved to be inconclusive. However, the K-S test 
appeared to have a higher power than the F^ test when the shape of the 
distribution was "near” normal such as the ones for 
(6 1 ,B 2 ) = (0,2.3), (0,2.8) and (.25,3.2). 
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TABLE 5 



n=21 



n-31 



n-51 



FOUTZ F TEST VS CHI SQUARED TEST 
n 

RANDOM STABILIZED STANDARD DISTRIBUTIONS 





0-0 


0.25 


0.50 


0.75 


1.00 


Sig 














Level 


a 














1.0 


102.0 


102.7 


101.6 


95.8 


96.9 


.10 


1.3 


96.1 


88.6 


69.1 


49.1 


11.7 


1.6 


72.0 


68.3 


59.5 


49.0 


29.2 




1.9 


54.3 


55.4 


55.4 


53.9 


55.8 




1.0 


120.2 


117.2 


121.6 


114.9 


110.3 


.05 


1.3 


97.9 


98.2 


71.2 


49.6 


- 1.4 


1.6 


71.2 


63.4 


58.6 


41.4 


20.1 




1.9 


49.8 


45.6 


51.1 


49.8 


44.2 




1.0 


256.0 


215.9 


238.0 


234.5 


233.8 


.01 


1.3 

1.6 


174.8 

118.5 


158.6 

130.3 


131.8 

107.9 


71.7 

62.1 


- 8.0 

21.7 




1.9 


81.6 


73.0 


99.0 


89.7 


78.7 




1.0 


37.9 


36.8 


38.1 


39.3 


37.0 


.10 


1.3 


35.9 


34.0 


28.4 


13.8 


-10.4 


1.6 


26.6 


26.2 


16.5 


11.1 


0.9 




1.9 


17.6 


19.5 


19.4 


13.4 


16.7 




1.0 


47.4 


44.6 


46.6 


50.2 


46.7 


.05 


1.3 


41.8 


38.9 


31.5 


11.7 


-18.6 


1.6 


28.8 


25.7 


13.0 


5.2 


- 7.6 




1.9 


12.0 


15.8 


13.3 


8.4 


14.4 




1.0 


74.7 


75.2 


73.8 


73.5 


70.2 


.01 


1.3 


55.8 


46.2 


31.2 


3.9 


-33.7 


1.6 


32.2 


30.7 


11.0 


- 6.0 


-21.3 




1.9 


14.5 


18.1 


10.1 


8.0 


1.0 




1.0 


11.4 


13.3 


11.1 


10.6 


11.4 


.10 


1.3 


9.6 


12.3 


6.1 


- 1.1 


-11.1 


1.6 


1.2 


- 0.5 


- 4.9 


- 7.3 


-13.9 




1.9 


-10.8 


- 8.9- 


-11.0 


-11.1 


-11.2 




1.0 


17.4 


18.5 


17.1 


16.3 


17.1 


.05 


1.3 


12.6 


12.2 


5.0 


—3.5 


-18.2 


1.6 


0.5 


- 3.5 


- 8.7 


-12.0 


-21.7 




1.9 


-14.3 


-12.9 


-15.1 


-16.3 


-15.9 




1.0 


30.5 


30.5 


31.2 


28.4 


25.8 


.01 


1.3 


14.1 


13.8 


3.1 


-12.2 


-35.2 




1.6 


- 5.4 


- 8.3 


-20.6 


-23.5 


-37.0 




1.9 


-24.2 


-24.2 


-24.1 


-24.2 


-26. L 
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TABLE 6 



FOUTZ F TEST VS KOLMOGOROV- SMIRNOV TEST RANDOM 
RANltf)M STABILIZED STANDARD DISTRIBUTIONS 



Sie 


B-* 0.0 


0.25 


0.50 


0.75 


1.00 




Level. 


a 
















1.0 


120.2 


123.8 


110.0 


120.4 


115.2 


* 


.10 


1.3 


102.8 


100.6 


76.0 


43.5 


0.9 






1.6 


87.7 


77.2 


61.3 


49.6 


22.7 






1.9 


65.5 


61.7 


62.9 


65.1 


57.8 






1.0 


197.5 


196.0 


186.7 


185.9 


192.1 




.05 


1.3 


170.8 


151.2 


111.3 


68.1 


- 2.8 


n=21 


1.6 


131.3 


110.4 


89.9 


67.0 


30.4 






1.9 


83.7 


90.5 


87.5 


85.6 


74.4 






1.0 


601.5 


465.5 


530.4 


522.7 


526.8 




.01 


1.3 


396.4 


338.1 


245.0 


121.4 


- 0.6 




1.6 


305.7 


243.9 


165.2 


117.7 


45.6 






1.9 


140.1 


147.5 


176.2 


176.5 


163.8 






1.0 


82.1 


80.6 


83.3 


84.7 


81.6 




.10 


1.3 


83.6 


77.6 


56.7 


25.5 


- 9.2 




1.6 


63.8 


59.5 


42.6 


29.4 


7.3 






1.9 


49.5 


55.6 


48.2 


44.8 


40.3 






1.0 


150.9 


149.8 


145.6 


152.7 


151.5 




.05 


1.3 


147.3 


129.2 


88.1 


40.5 


- 13.0 


n=31 


1.6 


107.8 


98.9 


65.0 


41.1 


11.4 






1.9 


77.2 


86.7 


69.5 


72.5 


68.2 






1.0 


428.5 


429.5 


447.0 


399.0 


471.3 




.01 


1.3 


386.1 


316.4 


190.6 


76.4 


- 16.2 




1.6 


267.3 


210.3 


136.9 


79.9 


24.1 






1.9 


190.8 


186.2 


153.0 


146.8 


136.8 






1.0 


38.5 


42.7 


38.9 


37.5 


37.7 




.10 


1.3 


49.3 


49.4 


29.6 


8.4 


- 10.3 




1.6 


40.0 


35.7 


24.0 


11.4 


- 6.7 






1.9 


24.4 


31.4 


26.1 


22.2 


19.5 






1.0 


77.6 


81.2 


78.0 


78.3 


76.0 




.05 


1.3 


94.4 


85.2 


53.2 


14.7 


- 16.4 


n=51 


1.6 

1.9 


79.4 

49.8 


61.4 

54.6 


37.9 

50.6 


17.0 

42.0 


- 8.1 

35.4 






1.0 


285.4 


287.0 


275.2 


270.9 


262.6 




.01 


1.3 


287.6 


243.0 


138.7 


37.7 


- 26.1 




1.6 


196.2 


173.3 


91.2 


41.0 


- 7.9 






1.9 


130.4 


126.3 


115.4 


111.8 


99.6 
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TABLE 7 



n=21 



n=31 



n=51 



FOUTZ TEST VS CHI SQUARED TEST 
MIXED NORMAL DISTRIBUTIONS 





Y+ 1.00 


0.30 


0.20 


0.10 


Sig 


N(0,a) 










Level 


- 












N(0,2) 


55.8 


69.3 


76.3 


40.7 


.10 


N(0,3) 


32.8 


86,8 


75.9 


44.7 




N(0 ,4) 


20.5 


101.4 


89.0 


75.7 




N(0,2) 


47.0 


46.4 


64.5 


23.7 


.05 


N(0 , 3) 


30.5 


71.3 


69.0 


34.7 




N(0,4) 


22.6 


91.1 


76.8 


48.1 




N(0,2) 


62.0 


262.1 


10 7.5 


88.6 


.01 


N(0,3) 


61.2 


109.8 


145.8 


86.3 




N(0,4) 


49.0 


118.3 


131.1 


125.0 




N(0,2) 


13.8 


46.9 


42.4 


22.0 


.10 


N(0, 3) 


4.3 


56.7 


50.6 


51.5 




N (0,4) 


1.8 


55.8 


58.7 


39.2 




N(0, 2) 


12.0 


38.3 


29.2 


28.6 


.05 


N(0 , 3) 


3.1 


53.2 


48.4 


53.2 




N(0,4) 


1.0 


57.4 


61.5 


44.4 




N(0 , 2) 


1.2 


62.1 


46.0 


52.1 


.01 


N(0 , 3) 


0.6 


47.3 


103.1 


104.7 




N(0,4) 


0.1 


84.0 


95.5 


61.5 




N(0,2) 


-13.1 


23.6 


26.1 


9.7 


.10 


N(0,3) 


- 5.9 


19.3 


28.4 


25.7 




N(0,4) 


- 1.4 


26.2 


29.5 


35.6 




N(0,2) 


-16.2 


28.9 


45.5 


17.4 


.05 


N(0 , 3) 


- 9.4 


29.2 


43.4 


38.0 




N(0, 4) 


- 2.5 


34.1 


45.4 


39.7 




N (0, 2) 


-26.9 


83.8 


54.8 


67.4 


.01 


N(0,3) 


-17.7 


41.6 


45.9 


67.9 




N(0,4) 


- 7.9 


24.5 


78.3 


51.4 
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TABLE 8 



FOUTZ F TEST VS KOLMOGOROV- SMIRNOV TEST 
n 

MIXED NORMAL DISTRIBUTIONS 





T* 


1.00 


0.30 


0.20 


0.10 




_ Sig 


N(0,a) 












Level 


N(0,2) ' 


60.5 


36.1 


31.8 


8.7 




.10 


N(0 , 3) 


52.6 


59.3 


42.9 


19.3 






N(0 , 4) 


38.8 


71.5 


54.0 


40.2 






N(0 , 2) 


87.4 


49.5 


36.4 


6.8 


n=21 


.05 


N(0, 3) 
N(0,4) 


91.8 

74.4 


75.4 

90.2 


49.8 

66.6 


16.9 

43.7 






N(0,2) 


156.3 


94.4 


43.1 


15.8 




.01 


N(0,3) 


223.6 


112.3 


71.0 


41.8 






N(0,4) 


198.1 


107.0 


101.2 


26.6 






N(0,2) 


42.6 


19.3 


27.8 


13.9 




.10 


N(0,3) 


31.3 


52.9 


38.8 


24.1 






N(0,4) 


18.5 


66.2 


52.9 


21.1 






N(0,2) 


69.3 


29.9 


27.8 


20.6 


n=31 


.05 


N(0,3) 
N(0 ,4) 


60.8 

38.7 


65.7 

101.4 


52.9 

61.9 


36.8 

36.2 






N(0,2) 


135.8 


52.9 


67.3 


55.3 




.01 


N(0,3) 


183.9 


107.5 


120.3 


66.0 






N(0,4) 


129.8 


163.9 


152.2 


72.1 






N(0,2) 


22.4 


26.8 


23.6 


9.7 




.10 


N(0, 3) 


9.6 


47.8 


39.4 


21.2 






N(0,4) 


3.1 


66.4 


51.1 


36.3 






N(0,2) 


45.7 


30.0 


39.9 


9.0 


n=51 


.05 


N(0, 3) 
N(0,4) 


24.5 

10.6 


74.0 

110.7 


52.2 

73.3 


29.9 

52.3 






N(0 , 2) 


117.6 


61.9 


81.1 


20.3 




.01 


N(0,3) 


96.0 


215.7 


83.3 


39.1 






N(0,4) 


53.0 

l 


170.9 


164.5 


86.7 
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TABLE 9 



FOUTZ F^ TEST VS CHI SQUARED TEST 
MIXED NORMAL DISTRIBUTIONS 



SIGNIFICANCE LEVEL 





MIXED NORMAL 


.10 


.05 


.01 


; 


0 . 7GxN (0.5, l)+0 . 30xN(0 .0,3) 


24.0 


- 0.7 


-20.8 


n=21 


0 . 80xN(0 . 5 , l)+0 . 20xN(0 .0 , 3) 


14.4 


- 9.2 


-27.2 




0 . 80xN (0 . 5 , 1) +0 . 20*N (0.0,4) 


17.2 


- 3.2 


-21.1 




0.70xN(0.5,l)+0.30XN(0.0,3) 


9.6 


- 6.8 


-28.6 


3 

II 

u> 

h-* 


0 . 80xN (0 . 5 , 1) +0 . 20xN (0.0,3) 


3.3 


-12.4 


-33.3 




0 . 80xN (0 . 5 , l)+0 . 20xN (0.0,4) 


3.7 


-10.7 


-30.2 




0 . 70xN (0 . 5 , l)+0 . 30xN (0.0,3) 


25.1 


-24.1 


-43.6 


n=51 


0 . 80xN (0 . 5 , l)+0 . 2GxN(0 .0,3) 


- 2.3 


-23.6 


-48.6 




0. 80xN(0 . 5 ,l)+0 . 20xN(0. 0,4) 


- 4.1 


-26.3 


-45.5 



i 
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TABLE 10 



FOUTZ F TEST VS KOLMOGOROV- SMIRNOV TEST 
n 

MIXED NORMAL DISTRIBUTIONS 



SIGNIFICANCE LEVEL 





MIXED NORMAL 


.10 


.05 


.01 




0. 70xN(0. 5 ,l)+0. 30xN (0.0,3) 


-17.0 


-28.9 


-33.7 


n=21 


0 . 80xN (0 . 5 , l)+0 . 20xN (0.0,3) 


- 27.3 


-38.2 


-40.7 




0 . 80xN (0 . 5 , l)+0 . 20xN (0.0,4) 


-25.5 


-32.5 


-35.1 




0. 70xN(0. 5 , l)+0. 30xN (0.0,3) 


-25.4 


-34.5 


-43.4 


n=31 


0. 80xN (0.5, l)+0. 20xN(0. 0, 3) 


-34.8 


-45.5 


-50.7 




0 . 80xN ( 0 . 5 , l)+0 . 20xN (0.0,4) 


-34.6 


-40.2 


-45.9 




0. 70xN(0.5,l)+0.30xN(0.0,3) 


- 28.7 


-47.6 


-56.7 


n=51 


0 . 80xN (0 . 5 , l)+0 . 20xN (0.0,3) 


-48.7 


- 57.4 


-64.9 




0 . 80 X N (0 . 5 , l)+0 . 20xN (0.0,4) 


-43.3 


-51.3 


-61.0 
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TABLE 11 









FOUTZ F TEST VS CHI 
n 


SQUARED TEST 










PEARSON DISTRIBUTIONS 














SIGNIFICANCE 


LEVEL 




'l 




S 2 


.10 


.05 


.01 




0.0 




2.30 


3.4 


-10.5 


-20.6 




0.0 




2.80 


21.1 


- 0.9 


-12.4 




0.01 




1.75 


12.8 


41.1 


57.8 




0.01 




1.90 


- 9.1 


- 2.4 


18.9 


n=21 


0.25 




3.20 


16.6 


- 0.4 


- 2.3 


0.50 




3.00 


17.5 


5.8 


- 9.3 




1.00 




3.40 


-12.5 


-16.2 


-15.8 




1.00 




3.60 


9.1 


-31.4 


-30.9 




1.00 




3.80 


16.1 


-17.6 


-24.8 




2.00 




5.50 


-16.1 


-14.9 


-20.7 




0.0 




2.30 


- 6.6 


-19.0 


-22.5 




0.0 




2.80 


10.7 


-12.1 


- 5.4 




0.01 




1.75 


- 5.3 


39.4 


83.5 




0.01 




1.90 


-18.7 


- 6.9 


26.9 


n=31 


0.25 




3.20 


0.7 


-11.8 


- 3.6 


0.50 




3.00 


- 1.1 


- 0.6 


-12.2 




1.00 




3.40 


-25.5 


-23.5 


-23.4 




1.00 




3.60 


- 5.7 


-41.4 


-41.7 




1.00 




3.80 


0.0 


-26.9 


-32.9 




2.00 




5.50 


-29.5 


-19.4 


-26.5 




0.0 




2.30 


30.4 


-39.0 


-19.7 




0.0 




2.80 


34.9 


-40.0 


-17.3 




0.01 




1.75 


6.6 


68.1 


132.9 




0.01 




1.90 


-11.8 


-15.8 


20.7 


n=51 


0.25 




3.20 


33.3 


-10.8 


-15.6 


0.50 




3.00 


2.5 


-19.9 


-20.2 




1.00 




3.40 


-37.4 


- 35.7 


-35.6 




1.00 




3.60 


-13.6 


-53.7 


- 56.7 




1.00 




3.80 


- 4.7 


-40.6 


-46.5 




2.00 




5.50 


-38.5 


-28.7 


-38.0 
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TABLE 12 



FOUTZ F TEST VS KOLMOGOROV- SMIRNOV TEST 
a 

PEARSON DISTRIBUTIONS 



SIGNIFICANCE LEVEL 





*1 


*2 


.10 


.05 


.01 




0.0 


2.30 


-22.8 


-25.0 


-26.2 




0.0 


2.80 


5.4 


-18.8 


- 7.2 




0.01 


1.75 


- 9.4 


20.6 


66.0 




0.01 


1.90 


-33.3 


-15.3 


16.5 


n=21 


0.25 


3.20 


- 8.1 


-14.8 


- 4.1 


0.50 


3.00 


9.9 


21.7 


34.9 




1.00 


3.40 


84.8 


94.4 


95.9 




1.00 


3.60 


52.9 


64.5 


76.7 




1.00 


3.80 


38.4 


44.0 


59.4 




2.00 


5.50 


62.4 


74.7 


70.4 




0.0 


2.30 


-28.9 


-36.8 


-36.9 




0.0 


2.80 


11.2 


-20.8 


-13.0 




0.01 


1.75 


-14.1 


27.6 


89.1 




0.01 


1.90 


-36.7 


-15.4 


19.4 


n=31 


0.25 


3.20 


-12.9 


-14.0 


- 9.3 


0.50 


3.00 


- 0.6 


25.0 


31.3 




1.00 


3.40 


100.0 


134.3 


131.7 




1.00 


3.60 


55.5 


79.4 


91.8 




1.00 


3.80 


39.9 


46.2 


65.3 




2.00 


5.50 


63.3 


91.0 


100.0 




0.0 


2.30 


- 9.1 


-48.6 


-40.0 




0.0 


2.80 


- 4.9 


-38.9 


-35.8 




0.01 


1.75 


-18.4 


62.0 


171.8 




0.01 


1.90 


-37.4 


-22.9 


26.6 


n-51 


0.25 


3.20 


- 8.1 


- 9.6 


-13.8 


0.50 


3.00 


0.0 


25.9 


47.0 




1.00 


3.40 


169.9 


225.8 


245.5 




1.00 


3.60 


64.9 


138.9 


154.4 




1.00 


3.80 


42.1 


69.0 


85.2 




2.00 


5.50 


103.0 


151.0 


164.7 
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5. CONCLUDING REMARKS 



The superior performance of the Foutz test in detecting certain types 
of deviations from the hypothesized distribution leads to several more prob- 
lems to be considered. Of primary importance is the generation of percentage 
points for the distribution of F^ for various values of n. The intractability 
of the problem of obtaining the exact distribution requires an empirical 
approach to finding a correction to the asymptotic approximation given by 
Foutz . 

Since the test is also applicable to p-variate distributions, an investi- 
gation of ways to obtain the statistically equivalent blocks, and then the 
probability content of them, at least for p = 2, is to be considered. The 
problems of obtaining these blocks and their contents becomes increasingly 
complicated in higher dimensions. 

The authors wish to acknowledge extensive use of the IBM 360/67 computer 
system at the Naval Postgraduate School, on which all computations were 
performed. 
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noon oooooooooooo 



APPENDIX 1 



this PROGRAM GENERATES DEVIATES tn random stabilized standard form 
BLGCKS GF 20, 30 » AND 50 DEVIATES ARE CONSIDERED BY GENERATING 50, 
THEN CONSIDERING THE FIRST 20, THE NEXT 30, AMD ALL 50. 

THE HYPOTHESIS THAT THEY CAME FROM A NORMAL DISTRIBUTION WITH 
MEAN AND VARIANCE AS GIVEN IN INPUT IS TESTED. 

THE TESTS PERFORMEC ARE THE CHI SQUARED, KOLMOGOROV- SMI RNOV , AND 
FOUTZ TEST AT THE CCNFIOENCE LEVELS OF 10, 5, AND 1 PERCENT. 

THE INPUT VALUES APE MEAN, VAR! ANCE, I SEED , AND THE NUMBER OF REPLI 
CATIONS. 

IMPLICIT R£AL*8 (A-Z) 

INTEGER** T, IER»N»I»NP1,N1»N2»N3*I SEED, NR , N STOP , IDF »IQ,ISEED1 

1 , NCEL ( 3 ) , NST ( 3 ) ,N SM P (3 i ,NCHI (3 ,3 ) ,NCHIN (3,3 ) , NKS (3, 3 ) , 

2 NKSN(2»3)»NFNN(2,3),NFN(3»3) ,NSR 
DIMENSION RD(55 J ,CHISQU(2) , 

1 KOLSMR ( 2 ) , C 0UTZL<2) ,FNTST(3,2) , ZA LF ( 3 ) , RN (55 I 
REAL** CELLS (50 l » CCM P ( 50 ) » Q» CH I , R ( 55) , PD I F ( 6) , QTST( 3 ) , RNN( 55) , 

1 STST( 3 » 3 ) 

DATA STST/. 26*73,. 21 756,. 165 5S,. 29*08, .2*170,. 188*1 ,.352*1, 

1 .28987,. 2260*/ 

DATA NCEL/8,10,10/,NST/l,21,l/,NSMP/20,30,50/,QTST/.l, .05, .01/ 

1 , CHIS QU / • CH I SQUA ' , * RED TE ST • /,KOLSMR/ • KOLMOG-S • , • MI R TEST*/ 

2 ,FOUTZL/' F0UT2 T E ST •/, ZALC1 , ZALC2/ . 2*306900, . 3678 79DQ / 

DATA WT1,MEAN1 , VAR1 , WT2 ,ME AN 2 , VAR2 /. 8D0 ,0. DO , 1 .DO , .2 DO ,0 . DO , 3 . DO/ 

3 , ZAL F/ 1. 281 55C0, 1.6**8500, 2.3263500/ 



THE NEXT * CARDS PUT 
FOUTZ DISTRIBUTION. 



IN THE EMPIRICAL CRITICAL VALUES FOR THE 



1 

2 



DIMENSION FN=mp<3,3) 

DATA FNEMP /. *271*00,. *190300,. *)816D0, 
.**865 CO , .*355300, .*211600, .*865900, 

• *6579 DO,. ***870 0/ 



100 



1* 

1 



17 

12 



110 

5 

3 



REAL** ALPHA, BETA ,BFRIME 
COMMON / NCR PRM/ ME AN » VARI » SVAR I 
EXTERNAL NORM , UN 1 F 
T = 1 

PIB2 = CAT AN(1.DC )*2.00 
PRINT 1 

READ 2, MEAN, VARI 
IF (MEAN+VARI .LE.O.CO)STOP 
FORMAT ( 2E5 .0 ) 

SVAR I * DSORT (VARI) 

WRITE ( * , 1* ) 

FORMAT ( 1H1 ) 

FORMAT ( '+INPUT MEAN AND VARIANCE FORMAT( 2E5. 0 ) * ) 

PRINT 17 

READ 12, ALPHA, BETA 

FORMAT ( ' ♦ INPUT ALPHA ANO BETA, FORMAT! 2F5.0) ' ) 

FORMAT ( 5F5 .0 ) 

PP INT 3 

READ 5 , ISEED»NB 

ISEED1 = I SEED 

BPRIME = BETA 

IF ( ALPHA .EC.O. )GC TC 110 

PHIZ = -PI B2*BETA*(1. - ABS ( 1 • - ALPHA )) /ALPHA 
BPRIME = DTAN(PIB2*(1. - ALPHA ) )*0TAN( PHI Z*ALPHA) 

CONTINUE 

DSEED = DFLOAT ( I SEED ) 

FORMAT (110,15) 

FORMAT ( • +INPUT SEED AND NUMBER OF REPLICATIONS, FORMAT! 110,15)* ) 

CO 205 IQ= 1 » 3 

DO 205 Nl*l,3 

NCHI(N1,IQ) = C 

NCHIN(NltlQ) = 0 

NKS(N1 , I C) = 0 

NKSN ( N1 , IQ ) = C 
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0 





c 

c 

c 

c 



THE FOLLOWING CARD PUTS IN THE EMPIRICAL VALUES FOR THE 



FOUTZ TEST 



FNTST (Nl , IQ ) = FNEMP ( Nl, IQ ) 



C 



205 NFNN ( N 1 , IQ ) = 0 
DO 300 1=1, NR 

CALL RA NSRV( 50 ,R N ,RD , CSEED , ALPHA , 8 PRIME ) 




N3 = NSMP(Nl) 




CALL GF IT (UN IF ,NC EL (Nl ) ,PNN , N3 , CELLS , COMP , CHI , IDF,Q, IER) 

DO 241 I C= 1 , 3 

I F ( Q • L T • GT ST ( 1 0 ) )NCHI(N1»IQ) =NCHI(N1,IQ) + 1 

241 CONTINUE 

CALL NKS1( UN IF ,RNN , N3 , PDI F , IS R ) 

CO 242 IC = 1,3 

IF(PQIF(1) .GT.STST(M,IQ))NKS(N1, IC) = NKS(N1,IQ) + 1 

242 CONTINUE 

CALL FGUTZ(UNIF,PNN»N3,Q) 

DO 243 I C=1 , 3 

IF(0.GT.FNTST(N1, 10) )NFN(N1, IQ) = NFN(N1,I0) + 1 

243 CONTINUE 
IDF = 0 

CALL GFIT(N0RM,NCEL(N1),R,N3, CELLS, COMP, CHI ,IDF,Q,IEP ) 

DO 244 IC = 1,3 

IF(0«LT.QTST(IQI )NCHIN(N1,IQ) = NCHIN(N1,I0) ♦ 1 

244 CONTINUE 

CALL NKS1(N0RM,R,N3,P0!F,I5R) 

DO 245 10=1,3 

IF(PDIF(1).GT.STST(M,IQ))NKSN(N1,IQ) = NKSN(N1,I0) ♦ 1 

245 CONTINUE 

CALL FOUTZ (NCR N,R,N2»Q) 

DO 246 I C=1 , 3 

IF(O.GT.FNTST( Nl ,IC) )NFNN(N1 , IC) = NFNN(N1,I0) + 1 

246 CONTINUE 
250 CONTINUE 
300 CONTINUE 

ISEED = CSEED 

WRITS (4, 7) MEAN,VAFI,NR,!SFED1, ISEEC 
WRITE! 4, 11 )ALP HA, BETA 

11 FORMAT ( // ' RANDOM STABLE STANDARDIZED FORM (RSSF), ALPHA = ',F5.2, 
1 *, BETA 5',F5.2//I 
WRITE (4, 6KHIS0U, ( ( NCHI (Nl , 1 0 ) , 

1 Nl = l,3) ,10=1,3) , ( ( NCHIN ( Nl , IC), 

2 Nl=l, 3) ,10=1,3) 

WRITE(4,8)K0LSMR,( (NKS(N1, IQ),N1=1,3) ,10=1,3) , 

1 ( <NK$N( M,IC) ,N1 = 1,3) ,10=1 ,3 ) 

WRITE! 4,8) FOUTZL, ( (NFN(N1, 10) ,N1 = 1 ,3) ,10=1,3 ) , (( NFNN (Nl, IQ) ,N1 = 1, 
1 ), 10=1,3) 

7 FORMAT (A/* THE VALLES OF MEAN (U) AND VARIANCE (V) ARE S2F12.5/// 
9 * THE RESULTS OF ',15,' P5PL ICATI CNS, STAR T ING WITH SEED', 

1 I12//33X, 'ENDING WITH NEXT SEEC',112//) 

8 FORMAT( ///20X,2A8//20X, '20 PIS', 

2 6X , ' 30 PTS ' ,6 X , '50 PTS V/10X, ' 10* ', 31 12/' CONTROL ', 2X , • 5?*, 
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END 



200 



SUBROUTINE RANSR V ( NO , RAN, R , D SEED , ALPHA » BPRIME ) 
IMPLICIT REAL*8 (A-M, 0-Z) 

DI MENSICN R ( 1 ) ,RAN <11 

REAL*4 RN< 100 > .ALPHA . BPRIME ,R STAB , Vs 

CALL GGUBS(DSEECt2*ND.RN) 

DO 200 NN= 1 » ND 

PN<NN+NDI ) 

R(NN) = RSTAB< ALPHA , BPRIME ,RN(NN) ,W> 

CONTINUE 

RETURN 

END 



SUBROUTINE NCRN(X.P) 

REAL*8 MEAN.VAPI .SVARI 
COMMON /NCR PR M/ ME AN tVARI .SVARI 
T = (X - MEAN I /S VAR I 
P = .5*£RFC(-T*.7G71068) 
RETURN 
END 



subroutine 

P = X 

RETURN 

END 



UNIF(X,P) 



FUNCTION RSTA6(ALPHA»8PRIMEtUtW) 

RSTAB IS A RANDOM STABLE STANCARCI ZED FORM 



(WHATEVER ) 



ARGUMENTS 
ALPHA = 
BPRIME = 
U = 

W 



50 



CHARACTERISTIC £XPCN£N T 
SKEWNESS IN REVISED PARAMETERIZATION 
UNIFORM VARIATE ON (0*1) 

EXPONENTIALLY DISTRIBUTED VARIATE 
DOUBLE PRECISION CA,CB 

DATA PI BY2/1. 57079633/, PI8Y4/ .785398 163/ *THR1/. 99/ 

EPS = 1. - ALPHA 
COMPUTE SOME TANGENTS 
PHIBY2 = PIBY2*<U - .5) 

A = PHIBV2*TAN2(PHI8Y2) 

B8 = TAN2(EPS*PH!EY2 ) 

B = EP S* PHI BY2*BE 

IF(EPS .GT.-.99ITAU = BPRIME /( TAN2 ( EP S*PI BY2 ) *P IBY2 ) 
IF(EPS.LE.-.99)TAL = BPRIME*P I BY2*5PS*( 1. -E°S ) *TAN2< < 1. 
1 PIBY2 ) 

COMPUTE SOME NECESSARY SUBEXPRESSIONS 
IF PHI NEAR PI BY 2* USE DOUBLE PRECISION. 

IF(A.GT.THR1 )G0 TO 50 
SINGLE PRECISION 
A2 = A**2 
A2P = 1. + A2 
A2 = 1. - A2 

B2 = 8**2 

B2P = 1. + B2 
B2 =1. - 82 
GO TO IOC 
DOUBLE PRECISION 
DA = CBLE(A)**2 
DB = DBLE( B ) **2 
A2 = l.CO - DA 
A2P = 1 • DO + DA 
B2 » 1. - OB 
D2P = 1. + DB 
COMPUTE COEFFICIENT 



- EPS)* 
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100 Z = A2P* ( B 2 ♦ 2.*PHIBY2*8B*TAU1/(W*A2*B2P) 

C COMPUTE THE EXPONENTIAL-TYPE EXPRESSION 

ALOGZ = ALOG(Z) 

0 = D2(EPS*AL0GZ/(1. - EPS ) )*( ALOGZ/( 1. - EPS)) 

C COMPUTE STABLE 

P.STAB = (1. + EP$*D)*2.*((A - e)*(l. + A*B) - PHI 6Y2 *T AU*BB*( B*A2 
1 - 2 .*A ) )/ (A2*E2P ) + TAU*0 
RETURN 
END 



100 



FUNCTION D2( Z ) 

EVALUATF ( EXP ( X I - 1)/X 

DOUBLE PRECISION F1,P2»G1,Q2,C3,PV,ZZ 

DATA P1»P2»Q1,Q2»Q3/. 84006685 2 5364 €3 23 903 ♦•23001114158996456902 ♦ 
. . 168013370507926648 C4, .1800 133704C7390023D3, 1 .DO/ 

THE APPROXIMATION 1801 FOP HART ET AL (1968, P213) 

IF ( ABS (Z).GT .0.1 ICC TO 100 

ZZ = Z*Z 

PV = PI 4 ZZ*P 2 

C2 = 2 . CO*PV/ ( 01 4 Z Z*(Q2 + ZZ*Q3) - Z*PV) 

RETURN 

D2 = (EXP(Z) - l.)/Z 

RETURN 

END 



FUNCTION TAN(XARG) 

C TANGENT FUNCTION 

LOGICAL NEG.INV 

DATA P0,PltP2,00,Ql ,02/ . 129 22103553 ,-.887662 3 77 El, .528644456E-1 , 
1 .1645 2 93 32 £3,- .451*2056 IE 2, 1./ 

C THE APPRCX IMATI CN 4283 FROM HART ET AL (1968, P. 251) 

DATA PIBY4 / . 785 39 8 1 6 3/ , PI B Y 2 / 1 . 57C79633 / , PI /3 . 1 4159265/ 

NEG = .FALSE. 

INV = .FALSE . 

X = XARG 
NEG = X.LT.Q. 

X = ABS( X) 

C PERFORM RANGE RECUCTION IF NECESSARY 

IF(X.LS.PIBY4)GC TC 50 
X = AMOD(X.PI) 

IF(X.LE.FIBY2)GQ TO 30 
NEG = .NCT.NEG 
X = PI - X 

30 IF(X.LE. FI BY 4) GC TC 50 

INV = .TRUE. 

X = PIBY2 - X 
50 X = X/PieY4 

C CONVERT TO RANGE CF RATIONAL 

XX = x *x 

TAN = X*(PO 4 XX*(F1 4 XX*P2I)/(C0 4 XX*<01 + XX*Q2») 

IF(NEG)TAN = -TAN 
IF(INV)TAN = 1 ./TAN 
RETURN 
END 



FUNCTION T AN2 ( XAR G ) 

C COMPUTE TAN(X)/X 

C FUNCTION DEFINED CNLY FCR AB S ( XARG ) . LE . P I BY 4 

C FOR OTHER ARGUMENT RETURNS TAN(X)/X, COMPUTED OIRECTLY 

DATA P0,P1 ,P2, CO, Cl, C2/. 129221 03 5E3, -.887662 377Elt .528644456E- 1 , 
1 .16452523 2c 3, -.45122056 15 2,1. / 

C THE APPRCX IMA T ICN 4283 FROM HART ET AL (1968, P. 251) 

DATA PIBY4,PIEY2,PI / .7853981 6 3 , 1 . 57079633 , 3 . 14 15926 5 / 

X « ABS ( XARG ) 

IF(X.GT.PIBY4)GC TC 200 
X = X/PIBY4 

C CONVERT TO RANGE CF RATIONAL 
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THIS PRCGRAM GENERATES MIXED NORMAL DEVIATES WT1*N ( MEAN 1 , VAR 1 ) + 
( 1 - WT1)*N!M£AN2 ,VAR2) . 

BLOCKS CF 20, 30, AND 50 DFEVIATES ARE CONSIDERED BY GENERATING 
50, THEN CONSIDERING THE FIRST 20, THE NEXT 30, AND ALL 50. 

THE HYPOTHESIS THAT THEY CAME FROM A NORMAL DISTRIBUTION WITH 
MEAN ANO VARIANCE AS GIVEN IN INPUT IS TESTED. 

THE TESTS PE RF CP. NEC ARE THE CHI SQUARED, KOLMOGOROV-SMIRNOV, AND 
FOUTZ TEST AT THE CCNFIOFNCF LEVELS OF 10, 5, AND I PERCENT. 

THE INPUT VALUES APE MEAN, VAR IANCE, ISEED, AND THE NUMBER OF REPLI 
CATIONS. 

IMPLICIT REAL* 8 !A-Z) 

INTEGER** T.IER.N.I ,NP1 ,N1 ,N2 ,N3,I SEEO.NR.NSTOP, IDF, IC, ISEED1 

1 , NCEL ( 3 ) » NST ! 3 ) »NSMP!3i*NCHI (3,3) ,NCHIN ! 3, 3) , NK S! 3 , 3 ) , 

2 NKSN< 3,3) ,NFNN( 3,3) »NFN (3 » 3 ) , NSR 
DIMENSION RD(55) ,CHISQU(2) , 

1 K0LSMR(2) ,FQUTZL(2) , FNTST ! 3 , 3 ) , ZALF ! 3 ) , RN t 55 ) 

REAL*4 CELLS< 50 ) ,CCM F < 50) , Q, CH I , R ! 55 ) , PO I F ! 6 ) , 

I STST ( 3 , 3 ) 

DATA STST/. 26473, .21756, .16959, .29*08, .2*170, 

1 . 28987, .226C*/ 

DATA NCEL/ 8, 10 , 10/ ,N ST/1, 2 1 , 1 / ,N$MP/20 , 30 , 50/ , QTST/. 1 ,.05,. 01/ 

1 , CHI SQL/' CHI SCUA* ,'RE0 TEST • /, KOLSMR/ • KOLMOG-S • , *MIR TEST'/ 

2 » FOUTZ L/* FOUTZ TE ST •/, ZALC1 , ZALC2/ .243069D0, .36787900/ 

DATA WT 1, MEAN1, V AR 1, WT2, MEAN 2, VAP 2/. 800,0.00,1. DO,. 200,0.00, 3. DO/ 

3 , ZALF/ 1. 281 55 CO, 1.6**8 500, 2. 3 263 5D0/ 



QTST ( 3 ) , RNN (55) 
188*1, .352*1, 



THE NEXT * CARDS PUT IN 
FCUTZ TEST. 



THE EMPIRICAL CRITICAL VALUES FOR THE 



C I MENSION FNEMP (3,3) 

OAT A FNEMP /. *271*00,.* 190 3D0,. 40 816DQ, 

1 .44865C0, .4355300, .4211600, .486 59D0 , 

2 .*657900, .4448700/ 

CQMMQN/NCRPRN/MEAN ,VARI ,SVARI 
EXTERNAL NORM , UN I F 
T = 1 

100 PRINT 1 

READ 2,MEAN,VARI 
IFIMEAN+VARI.LE.O.CO )STOP 

2 FORMAT ( 2E 5 .0 ) 

SVARI = CS CRT ( VAR I ) 

WRITE (4, 14) 

14 FORMAT ( 1H1 ) 

1 FORMAT ! ' +1 NPUT MEAN AND VARIANCE, FORMAT ( 2F5.0 )• ) 

PRINT 17 

READ 12»WT1,MEAN1 , VAR1, MEAN2 , VAR 2 
WT2 = 1 CO ~ WT1 

17 FORMAT! »+ INPUT WT1 , M5AN1 ♦ VA R1 ,MFAN2 ,VAR2 , POP.MAT <5 F5 .0 ) • ) 

12 FORMAT! 5 F5.0) 

PRINT 3 

READ 5, ISEED, NR 
ISEED1 = ISEED 
DSEED = CFLOAT(ISEEC) 

5 FORMAT! 110,15) 

3 FORMAT! '+INPUT SEEC AND NUMBER OF REPLICATIONS, FORMAT! 110,1 5) ' ) 
CO 205 10*1,3 

DO 205 N 1= 1, 3 
NCHI ( N1 , I C ) = 0 
NCHIN!N1,I0) = 0 
NKS1N1, IQ) = 0 
NKSN ! N1 , IQ ) = 0 
NFN! N1 , IQ) = 0 

FNTST ! N1 , I C ) = ZALC 1 *ZALF ( IQ ) /OSORT! DFLOAT! NSMP! N1 )+ 1 ) ) + ZALC2 
FNTST ! N1 ,1 C) = FNEMP !M,IQ) 

205 NFNN! N 1 , IQ ) = 0 
DO 300 1=1, NR 
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CALL M I XNRM( 50,RN,RD ,DSSFD,WT1 ,MEAN1 , VAR1 , ME AN2 , VAR2 ) 

CO 250 M = l»3 
N2 = NST(Nl) 

N3 = NSMP(Nl) 

CALL VSRTAD(RD(N2),N3) 

CALL VSRTAD(RN(N2),N3) 

CO 210 NSR=1»N 3 
RNN(NSR) = RN(NSR -» N2-1) 

210 R( NSR l = RD( NSR+N2-1 ) 

IDF = 0 

CALL GFITCUNIF ,NCEL ( N1 ) ,RNN, N3 ,CELLS , COMP , CHI , IDF,Q, I5R) 
CO 241 10=1,3 

IF (Q.LT .CTST ( IC ) INCH I (Nl, IQ) ..= NCHI(N1»I0) + 1 

241 CONTINUE 

CALL NKS1(UNIF,RNN,N3,PDIF,!ER) 

DO 242 IC * 1.3 

IF(P0IF(1) .GT.STST(M,IQ) )NKS(M,IC) = NKS(N1,IG) + 1 

242 CONTINUE 

CALL FGUTZ(UNIF ,PNN,N3,Q) 



DO 243 10=1,3 

- - f$T(M, 
E 

IDF = 0 



IF(O.GT.FN 
243 CONTINUE 



IC) >NFN(N1, IQ) = NFN(M1,IQ) + 1 



CALL GF IT ( NORM, NCEL ( N1 ) ,R , N3 , CELLS ,COMP ,CHI , I DF* Q, IER ) 

DO 244 IC = 1,3 

IF(O.LT.QTST<IQ) INCH INC N1 , 1 Q ) = NCHIN(N1,IQ) + 1 

244 CONTINUE 

CALL NK$1(N0RM,R ,N3 .PDIF.ISR) 

CO 245 10=1,3 

IF(PDIF(1).GT.STST(M,IQ))NKSN(N1,I0) = MKSN(N1,I0) + 1 

245 CONTINUE 

CALL FOUTZ (NCR M, R , N3 , Q ) 

DO 246 10=1,3 

IF ( 0 ,GT * P NTST( N1 , 1 C ) )NFNN(N1,IC) = NFNN(N1,I0) + 1 

246 CONTINUE 
250 CONTINUE 
300 CONTINUE 

ISEED = CSEED 

WRITE ( 4, 7) ME AN, VAR I , NR, I SEED 1 ,ISEEC 
WRITE (4, 11 )WT1,MEAN1,VAR1,WT2,MEAN2,VAR2 
11 FORMAT ( //' MIXED NORMAL (MIXNCPM) = • , F5 .2 , • *N ( • , F3 . 1 , • , • , F 3 . 1 , 

1 ') ♦ ' , F5. 2, '*N(',F3.1»','»F3.1,' )'//) 

WPITE(4, 8)CHISCU, ( (N CHI (Ml, IQ) . 

1 Nl=l,3) ,10 = 1,3) ,((NCHIN(N1,IC), 

2 Nl=l, 3 ) , 10=1, 3 ) 

WRITE(4,8)K0LSMR,( (NKS(N1,!0) ,N1=1,3), 10=1,3), 

1 ( (NKSN(M,IO) ,N1 = 1 ,3) ,10=1,3) 

WRITE(4,8IF0UTZL. ( (NFN(N1, IQ) ,N1=1,3) ,10=1,3) , ( ( NFMN ( N1 ,1 0 ) ,M1=1 
1 ) , 10 = 1 ,3) 

7 FORMAT ( // ' THE VALUES OF MEAN (U) AND VARIANCE (V) ARE* ,2F12.5// 



9 • THE RESULTS OF ',15,' REPLICATIONS, STARTING WITH SEED', 



1 I12//33X, 'ENDING WITH NEXT SEED', 112//) 

8 FORMAT (///20X, 2A6//2CX, '20 PTS', 

2 6X , ' 30 PTS' ,6 X , • 50 PTS • //10X , • 10 % • , 31 12/ • CONTROL ', 2X, • 5T', 

3 31 12/10X, ' 1*',2I12///' MIXNCPM' ,2X, '10*' ,3112/' VS',5X,» 5*' 

4 3112/' N (U ,V I ' , 2X , • IX', 3112) 

WRITE ( 4,9) (( FNTSTCNl ,10) ,N1 =1,3), IC=1,3 ) 

9 FORMAT ( / '0 THIS WAS BASEC ON THE FOLLLOWING TEST VALUES', 

1 / (13X , 3 F12.5 ) ) 

GO TO ICO 
END 



SUBROUTINE M IXNRM (NO , RAN,R ,DSEED ,WT1 ,MEAM1 , VAR 1 ,MEAN2 , VAR 2 ) 
IMPLICIT REAL*8 (A-M,0-Z> 

DIMENSION R( 1) ,RAN( 1) 

REAL*4 RN(IOO) , RNN » RNNN 
SVAR1 = DSQRT(VAFl) 

. SVAR2 = DSQRT ( VAR 2 ) 

WT2 = l.CO - WT1 
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XX = X* X 

TAN2 = (P0 + XX*(P1 ♦ XX*P 2) ) /(PIBY4*(QQ + XXMQ1 + XX*02))» 
RETURN 

TAN2 = TAN(XARG) / XA PG 

RETURN 

END 



200 

300 

400 



SUBROUTINE FCUTZ (FCCF,X T »NXT»FN) 

THIS SUBROUTINE GENERATES THE STATISTIC FOR THE FCUTZ 
INPUT VARIABLES ARE 



FN TEST, 



PCDF - THE CUMULATIVE DISTRIBUTION FUNCTION AGAINST WHICH THE 

DEVIATES ARE BEING TESTED. CALLING SEQUENCE MUST BE OF 
THE FCRM • CALL PCDF(X,P) ', WHERE X IS AN INPUT VALUE, 
AND THE VALUE OF THE CUMULATIVE DISTRIBUTION FUNCTION 
IS RETURNED IN P. 

P MUST EE BETWEEN 0 AND 1. 

XT - THE ARRAY OF DEVIATES, IN INCREASING ORDER. 

NXT - THE NUMBER OF DEVIATES IN T H5 ARRAY XT < * N - 1 > 

THE RETURNED VALLE IS FN, THE VALUE OF T HE STATISTIC. 



■NT LY LIMITED 
!TD. 



NXT IS PRES_ 

THE VARIABLE X 
CIMENS ION XT ( 1 ) 

REAL*8 XTD ( 51 ) ,P N , FN D 
N = NXT ♦ 1 
DO 200 1 = 1, NXT 
K = N - I 

CALL PCDF(XT(K)»P) 

XTD ( K+l ) = P 

RN = l.DC/N 

XTD(l) = RN - XTD ( 2 ) 

DO 300 I =2, NXT 

XTD(I) = RN - XTC ( I +1 ) + XTD ( I ) 
XTD(N) = RN - l.CC + XTD( N ) 

FND = 0. 

CO 400 I = 1 ,N 

FNO = FNC ♦ CMAX1(XTC(I),0.00) 

FN = FND 

RETURN 

END 



TO A MAXIMUM OF 50 BY THE DIMENSION OF 
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150 

200 



CALL GGUBS(DSEED,2*ND,RN) 

CO 200 NN=1 , NO 
RAN(NN) = PN(NN) 

CALL MDNRI$(RN(NN),RNN,IER) 
IF(RN(NN4ND).GT.WT2 )G0 t 0 150 
R(NN) = PNN*SVAR2 + MEAN2 
GO TO 200 

R ( NN ) = RNN+SVAR1 + MEAN1 

CONTINUE 

RETURN 

END 



SUBROUTINE NCRR(XtP) 

PEAL*8 MSAN.VARI ,SVARI 
COMMON/ NOR PRM/ME AN »VARI,SVAPT 
T = (X - MEANI/SVAPI 
P = .5*£RFC<-T*. 1071068) 
RETURN 
END 



SUBROUTINE UNI F < X , P I 

P * X 

RETURN 

END 



200 

300 

400 



TS 



EST. 



SUBROUTINE FOUTZ ( PCOF,XT,NXT , FN ) 

THIS SUBROUTINE GENERATES THE STATISTIC FOR THE FCUTZ FN 
INPUT VAPIABLES APE: 

PCD F - THE CUMULATIVE DISTRIBUTION FUNCTION AGAINST WHICH THE 
DEVIATES ARE BEING TESTED. CALLING SEQUENCE MUS T RE 
THE FORM • CALL PCCF(X,P| », WHERE X IS AM INPUT VALU 
AND ThE VALUE OF THE CUMULATIVE DISTRIBUTION FUNC T ION 
IS RETURNED IN P. 

P MUST eE EETWE C N 0 AND 1 . 

XT - THE ARRAY OF DEVIATES, IN INCREASING ORDER. 

NXT - THE NUMBER OF DEVIATES IN THE ARRAY XT ( * N - 1) 



THE RETURNED VALLE IS FN, T H£ VALUE OF THE STATISTIC. 



NXT IS PRESENTLY LIMITED TO A MAXIMUM OF 50 BY THE DIMENSION OF 
THE VARIABLE XTD . 

DIMENSION XT (1 ) 

REAL*8 XTD ( 51 ) ,RN ,FNO 
N = NXT 4 1 
DO 200 1=1, NXT 
K = N - I 

CALL PCCF(XT(K),P) 

XTD(K+1> = P 

RN = l.DO/N 

XT D ( 1 ) = RN - XTC ( 2 ) 

DO 300 I =2, NXT 

XTD(I) = RN - XTO ( I +1 ) 4 XTD ( I ) 

XTD(N) = RN - l.CO 4 XTC(N) 

FND = 0. 

CO 400 1=1, N 

FND = FNC 4 CMAX1 <XTD(I),Q.DO) 

FN = FND 

RETURN 

END 
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THIS PROGRAM GENERATE S P? ARSON TYPE I OR II RANDOM DEVIATES. 

THE PARAMETERS ARE CALCULATED IN TERMS OF 81 AND B2 IN 
SUBROUTINE PRM. SUBROUTINE ACINT1 IS US K D tq CALCULATE THE 
CDF TABLE. GGUBS IS USED TO GENERATE RANDOM (O.ll DEVIATES 
(UNIFORMLY DISTRIBUTED, AMD THEN SUBROUTINE RANDP1 DOES 
AND INVERSE CALCLLATICN TO OBTAIN THE RANDOM 
PEARSON CEVIATE.) 

A SET OF 50 RANDOM CEVIA T ES ARE GENERATED, THEN THE CHI 
SQUARED TEST, THE KCLMOLGOROV-SMI RNOV TEST, AMD THE FOUTZ <=N 
TEST ARE APPLIED. THE UNIFORM DEVIATES ARE TESTED AGAINST 
THE HYPOTHESIS THEY CAME FROM A UNIFORM DISTRIBUTION AS A CON T ROL 
THEN THE PEARSON OEVIATES ARE TESTED AGAINST THE HYPOTHESIS 
THAT THFY CAME FPCM A NORMAL DISTRIBUTION. 

THIS TEST IS REFLICATEO A NUMBER OF TIMES. 

THE TESTS ARE APPLIED TO THE FIRST 20, t H c NEXT 30, AND ALL 50. 



THE SET CF DEVIATES IS ALSO TESTED AGAINST THE 
THIS IS REPEATED FOR THE SAME SETS. 



NORMAL DISTRIBUTION 



THE INPUT VALUES ARE B1.B2, I $E C 0, AND THE NUMBER OF REPLICATIONS 
IMPLICIT REAL*8 (A-Z) 

INTEGERS T,IER,N, I , NP1,N 1 , N 2 ,N3 , 1 SEED , NR ,NSTOP , I DF , I Q , I SEED! 

1 ♦ NCEL ( 3 ) , NST ( 3 ) , N$MP (3 ) , NCH I ( 3 ,3 ) , NCHIN ( 3, 3 > , NKS ( 3, 3 > , 

2 NKSN( 3 » 3 1 »NFNN( 3 ,3) ,NFN(3,3) , NAFNN ( 3 ,3 ) , NA'N ( 3 , 3 ) , NS R 

DI MENS ION X( 2001 ) ,COF( 2001 ) , A ( 2001 ) , R ( 2001 » , RD ( 55) ,CHISQU(2) , 

1 KOLSMR ( 2 ) » P 0UTZL(2) ,ASFL(2) , f=NTST (3 ,3 ) , ZALF ( 3 ) , RN( 5 5 ) 

COMMON /PPARM/C0,C1,C2, Al, A2,BCC,BC1,M1 ,M2,K!NV,XL,XR,MEAN,T 
COMMON /CDF/X,CCF,A,B,NP1 

REAL* A CELLS (50) ,CCM P ( 50) , Q , CHI , P ( 55 ) , PD I F (6 ) , QTST ( 3 ) , RNN ( 55 ) , 

1 STST (3,2) 

DATA STST/. 26473,. 21 756,. 16959,. 29408, .24170, . 18841, . 35241, 

1 .28987, .22604/ 

CATA NCEL/ 8, 10 , 10/, NST/1, 21 , 1 / ,NSMP/ 20, 30 , 50/ , QTST/. 1 ,. 05 , . 01 / 

1 ,CHISCU/'CHI SQUA' , ‘RED TEST*/. KCLSMR/ »KOLMOG-S • , 'MIR TEST'/ 

2 , FOUTZ L / ' EMP FOLT'.'Z TEST* /, ZA LC 1 , ZALC2/ .243 06 9C0 , .367879 00/ 

3 , ZALF/1.28155C0, 1.6 448 5D0, 2. 3 26 3 SCO/ 

DIMENSION FNEMP (3 , 3 ) 

DATA FNEMP/. 42714D0, .4190300, .4081600 , 

2 .4436500, .4355300, .4211600 , 

3 .4865900, .4657900, .4448700/ 

4 , AS FL/ ' AS YM FOU'.'TZ TEST’/ 

EXTERNAL NORM , UN I F 

T = 1 

100 PRINT 1 

READ 2,B1,B2 
IF(B1+B2.LE.0.D0) STCP 

2 FORMAT (2E5.0) 

WRITE (4,14) 

14 FORMAT ( 1 H 1 1 

1 F0RMAT(»4 B1 AND B2, FORMAT ( 2= 5.0 )• ) 

CALL PRM(B1,82,1.C0,IER) 

IF ( IER . NE . 0 ) GO TG SCO 
A(l) = 0.00 
CDF(l) = 0. 

CALL INTSZ(N.X) 

NP1 = N ♦ 1 
PRINT 3 

READ 5, ISSED.NR 
ISEED1 = I SEED 
DSEED = CFLOATd SEED) 

5 FORMAT (110, 15) 

3 FORMAT ( * +1 NPUT SEED AND NUMBER OF REPLICATIONS, FORMAT( 110, 15 ) • ) 
CO 205 I Q= 1, 3 

DO 205 Nl=l , 3 
NCHI (Ml ,IQ) = 0 
NCHIN(NI.IQ) = 0 
NKS (Ml, I C) = 0 
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205 



210 



241 



242 



243 



244 



245 



246 

250 

300 



+ ZALC2 



NKSN(NltlO) = 0 
NFN ( N1 1 10=0 

FNTST< N1 *10) = ZALC1*ZALF(10) /DS0RT(DFL0AT(NSMP(N1)+1) ) 

NAFNN< N1,IQ) = 0 
NAFN(N1,IQ) = 0 
NFNN<N1,IQ) = 0 
DO 300 I =1 »NR 

CALL RANPD1(50»RN,RC,0SSFD) 

DO 250 Nl=l» 3 
N2 = NST(Nl) 

N3 = NSMP(Nl) 

CALL VSPTAD<RD(N2) ,N2) 

CALL VSRTAD(RN(N2 ),N2) 

DO 210 NSR=1 »N3 
RNN(NSR)= RN ( NSR ♦ N2-1) 

R( NSR I = RD(NSR+N2-1> 

IDF = 0 

CALL GFIT(UNIF,NCEL(N1),RNN,N3,CELLS,C0MP,CHI , IDF t Q t I ER) 

CC 241 10=1,3 

IF(0,L7.0TST<IQ) INCHI(N1,IQ) = NCHI(N1,IQ) + 1 
CONTINUE 

CALL NKSKUNIF , RAN, N3, PDIF, IER) 

DO 242 10 = 1,2 

IF(P0IF(1I .GT.STST(Nl,IO)NKS(Nl, IQ) = NKS(N1,IQ) + 1 
CONTINUE 

CALL F0UTZ(UNIF»RNN»N3»Q) 

DO 243 I C=1 r 3 

IF<Q.GT.FNEMP<N1,IC) )NFN(N1,IQ> = NFN(N1,I0> + 1 
I F (0, GT • FNTST ( N1 , 10 ) ) NA p N ( N1 , IQ ) = NAFN< N1 » I Q) + 1 
CONTINUE 
IDF = 0 

CALL GFIT(N0RM,NCEL(N1),R,N3, CELLS, COMP, CHI , I DP, Q, IER ) 

DO 244 IQ = 1,3 

IFIO.LT .OTST(IQ) )NCHIN(N1,!Q) = NCHIN( N1 ,10) + 1 
CONTINUE 

CALL NKS1(N0RM,R,N3, PDIF, IER) 

DO 245 10=1,3 

IF(PDIF(1).GT.STST(M,IQ))NKSN(N1, IQ) = NKSN(N1,IQI ♦ 1 
CONTINUE 

CALL F0UTZ(N0RM,R,N2,Q) 

DO 246 10=1,3 

IF(0.GT.FNEMP(N1,IC))NFNN(N1,!C) = NFNN(N1,I0) + 1 
IF(0.GT*FNTST(N1,IC) )NA P NN ( Ml , 10 ) = NAFNN( N1 , IQ) + 1 
CONTINUE 
CONTINUE 
CONTINUE 
ISEED = DSEED 
WRITE (4 ♦ 7 ) B1 , B2, A 1, A 2, Ml, M2, NR, I SEED 1,1 SEED 
WRITE(4,8ICHIS0U,< (NCHI (N1 ,10) , 

1 N 1=1, 3 ) ,10=1, 3 ) , ( (NCHIN(NI.IQ), 

2 N 1 = 1 , 3 ) ,10=1,3) 

WRITE(4,£)KCLSf'R,( (NKS(N1,I0) , Nl=l ,3 ) , 10=1 ,3 ) , 

1 ( <NKSN(N1,IQ) ,N1=1,3),IQ=1,3) 

WRITE (4,81 FOUTZL, ( (NFN(N1, IQ) ,N1=1,3), 10=1,3), ( ( NFNN ( N1,IQ ) ,N1=] 
1 ), 10=1,3) 

WR ITE (4 , 8 ) AS FL , ( ( NAFN <N1, 1 0 ) , Nl= 1, 2) , IQ = 1, 3) , 

1 ( (NAFNN(N1,IC) , N1 =1 ,3 ) , 1 0=1 ,3 ) 

7 FORMAT ( / / • THE VALUES OF B1 AND B2 ARE • ,2F12 .5 /// 

A • VALUES 0= Al, A2, Ml, AND M2 ARE* //4F12. 6/// 
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1 



• THE RESULTS CF 
I12//33X, ‘ENDING 



‘ , 1 5 , ‘ REPLICATICMS, STARTING WITH SEED' 



WITH N9 XT SEED* ,112//) 



8 FORMAT <///20X,2A8//2'3X,‘20 PTS*. 



,6X, ' 5C PTS'//10X,'10X' ,3112/' 

1** ,2112///' PEARSON ' , 2X , • 10* ',31 12/' 
IX' ,3112) 



2 6X,'30 PTS 

3 3I12/10X, • 

4 3112/' 

GO TO 100 

900 PRINT 4,N,IFR,B1,E2,ERR0R 

4 FORMAT ( ' N »I5R, 61 »B2 ' ,2 110 , 1 P3 C15 .6 ) 
STOP 
END 



VS' ,5X 



5* 
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SUBROUTINE RAN PCI (NIC ♦ RAN.R t DSEED ) 
IMPLICIT REAL* 6 (A-M.C-Z) 

CIMENS ION R ( 1 ) »R AN < 1 ) 

REAL*4 RMlOOItRNN 
CALL GGUBS?DSEED ,NC,RN) 

DO 200 NN=1»ND 
RAN< NN I = RN(NN) 

CALL CDF INV( RAN( NN ) t RNN) 

R(NN) = PNN 
CONTINUE 
RETURN 
END 



SUBROUTINE NORMX.F) 

P = .5*ERPC(-X*.7C71C68) 

RETURN 

END 



SUBROUTINE UNI F(X,P) 

P = X 

RETURN 

END 



SUBROUTINE PRM ( B1 t E2 tU2 » IER ) 

IMPLICIT REAL* 8 (A-Z) 

COMMON /PPARM/C0.C1.C2»A1, A 2 » BIGCO »BI GC1 , Ml , M2 t K IN V, XL ♦ XR » MEAN » T 
INT£GSR*4 T»IER 
IER = 0 

I F (T.EO • 2 * 82 = ( 6 . + 3.*Bl)/2. 

IF ( T.E C . 5 ) 31 = 4.* <2. *82 - 6 . ) * (4 . *B2-3 . I / < ( 82+3 . )**2 + 

1 12.*<4.*B2-2. 1 1 
DEN = 10.*B2 - 12.*ei - 18. 

CO = ( 4 • *B2 - 3.*ei)*U2/DEN 
Cl = DS0PT(U2*ei )*<B2 + 3.)/DEN 
C2 = (2.*B2 - 3.*E1 - 6 • ) / DEN 
A1 = 0. 

A2 = 0. 

BIGCO = 0. 

BIGC1 = C. 

Ml = 0. 

M2 = 0. 

GO TO (100,200,30C,400,500,600,7CO)tT 
100 CALL RCCTS (CO » Cl* C2. Al, A2, IER) 

IF ( IER. NE.O) RETURN 
Ml = (Cl ♦ A1)/C2/(A2 - Al) 



M2 = -(Cl + A2)/C2/(A2 - Al ) 
KINV = <A2 - A 1 I **( M1+ M2 + 



1. )*DGAMMA<M1*1. )*D GAMMA (M2 +1 . ) 

1 /DGAMMA (Ml + M2 * 2.) 

XL = Al 
XP = A2 
GO TO 750 
200 GO TO 100 
300 Ml = (C0/C1 - Cll/Cl 
C2 = 0. 

KINV = DGA MMA( Ml + 1 • ) A DEXP (C0/C1**2 )*DABS ( Cl ) **( 2 . *M1 + 1 . ) 

XL = -C0/C1 
XR = 1.C50 
GO TO 75C 

BIGCO = CO - Cl**2/C2/4. 

BIGC1 = C1/C2/2. 

Ml = (Cl - BIGC1 )/0$QRT(2. DO) /BIGCO 
M2 = DSQRT (BIGC0/C2 ) 

KINV = 1.00 
XL = -1.D50 
XR = 1.050 
GO TO 75C 

500 BIGC1 = C1/2./C2 



400 
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600 



Ml = (Cl - BIGCD/C2 
IF(C2.GT.1.0C)G0 70 800 

KINV = CABStMl )**(1./C2 - 1 . > /DGAMMA ( 1. /C 2-1. ) 

XL = -BTC-C1 
XR = 1 .050 
GO TO 750 

CALL R0CTS(C0*Cl»C2tAltA2»IER) 

Ml = -(Cl ♦ A 1 ) /C 2 /( A2 - Al) 

M2 * (Cl ♦ A2)/C2/(A2 - Al I 

IF(M2.GE.-1.D0.0R.M1 + M2 . G5 .0 .00 » GO TO 800 

KINV = ( A2 - A1)**(M1 + M2 + 1. ) *D GAMMA ( M 2+1. ) *DGAMMA ( -Ml -M2-1 • ) 

1 /DGAMMA (-Ml) 

XL = A 2 
XR. = 1.050 
GO TO 750 

700 KINV=C0**(-.5D0/C2)*DSQRT(C0/C2)*DGAMMA(.5D0)*DGAMMA< .5D0/C2-.5C 
1 /D GAMMA (.5D0/C2 ) 

Cl = 0. 

XL = -1.050 
XP = l C50 

MEAN = * ( 4. D0*B 2 - 3 .D0*B1 ) /3. DO/ ( B2 - B1 - 1 . DO ) *U2*PCF(0 .00 ) 
RETURN 
IER = 1 
RETURN 
END 



750 

800 



500 



600 



SUBROUTINE ROOTS ( CO, C1,C2, A 1 , A 2, IER) 
IMPLICIT REAL* 8 (A-Z) 

INTEGER*4 IER 
IER * 0 

IF(C2.EC.0.D0)GC TO 500 

DIS = C 1**2 - 4. CC*C0*C2 

IF(DIS.LT.0.D0»GC TC 600 

SDIS = 0 SORT ( D I S I 

DNUM = -Cl - SDIS 

IF(C1.LT.0.D0)CNUM = -Cl + SCIS 

X2 = DNUM/2.DO/C2 

XI = C0/C2/X2 

Al = DMIM(X1,X2) 

A2 = DMAX1 ( Xl» X2 ) 

RETURN 

A2 = 1.D75 

Al = -C0/C1 

RETURN 

IER = -1 

RETURN 

END 



FUNCTION PDF (X ) 



THIS FUNCTION EVALUATESTHE PEARSON DISTRIBU T ION 
FOR A GIVEN X, THE PARAMETERS HAVING PREVIOUSLY 8EEN 
CALCULATED IN SLBRCLTINE PRM. 



100 



IMPLICIT REAL*8 (A-2) 

INTEGER*^ ITY 

COMMON/ PPARM/ CO t Cl t C2» A 1 , A 2 » B IGCO t BI GC 1 , M 1 »M2 »K I NVt XI , X2 t MEA M « I 
GO TO (lQOt200»300»400t500»600t70Q)t!TY 
IF( X.LE.A1.0R.X.GE.A2)GO TO 140 
PDF = (X - A1I**M1*(A2 - X )**M2/K INV 
RETURN 

140 IF(X.LE.A1.AND.M1.LT.0.D0»G0 TO 150 
IF(X.GE.A2.AND.M2.LT.0.C0)G0 TO 15C 
PDF = 1.0-25 
RETURN 

150 PDF = 1.C25 
RETURN 

200 C-0 TO IOC 



44 



300 PDF = (CO + C1*X )**I*1*DEXP(-X/C1 )/KlNV 
RETURN 

400 PDF = (BIGCO + C 2* ( X+BIGC 1 ) **2 ) ** ( -1. /C 2) * 

1 DEXP ( -M*DATA N ( (X + P IGC1 ) / M2 ) ) / K INV 
RFTURN 

500 PDF = (X + BIGC1 )**(-l./C2 )*DSXP(Ml/( X ♦ BIGCDI/KINV 
RETURN 

600 PDF = (X - A 1 ) 1* ( X - A2)**M2/KINV 
RETURN 

700 PDF = (CC ♦ C2*X**2)**(-.5D0/C2)/KINV 
RETURN 
END 



150 

160 



170 



300 



301 



302 



304 



308 



309 

310 



400 



SUBROUTINE INTSZ(N.X) 

IMPLICIT REAl*8 (A-Z) 

COMMON /PPARM/C0»C1,C2, Al ? A2t 8C0tBCl l Ml,M2tKINVt XL,XR,M=AN t T 
COMMON /CCF/XDUMC 2001) , CDF ( 2001 ), A ( 2001 ), B ( 2001 ), NP1 
Cl MENS ICN X(l) 

INTSGER*4 T.N, NN , I NT , IEND » NP1 , IER 
X(l) = A 1 
SCI) = PCF(Al) 

CDF ( 1 ) = O.DO 
INT = 1 
NN » 1 

DXM = (A2 - A1J/1.C7 
CX = DXM*1 *03 
GO TO 160 

DX = D X* ( 1 .DO +DSCRT( .475D-4/-REST ) ) /2. 

INT = 1 

XCNN+l ) = X(NN) 4 CX 
IEND = 1 

IF(X(NN + 1) .LT.A2 )G0 TO 170 

X( NN+l ) = A2 

DX = X ( NN+l) - X ( N N ) 

IEND = 2 

B( NN4-1 ) = PDF( X( NN4-1 ) ) 

ESTINT = ( B ( NN 4 1 ) 4 B (NN ) ) *D X/ 2. DO 

IF(ESTINT.LT.l .C-8 )GQ TO 301 

DC DF = ADINTK X(NN) ,X(NN4-1) t FRRCP.»ISR) 

IF ( IER. GT. 100) GO TO 900 
IF(DCDF.GT.,025D0)GC TO 309 

EREST = (l.D0/B(NN4l) - l.DO/8 (NN) )*DCDF/8 .DO 
A( NN) = EREST 
EREST = CABS(EREST) 

GO TO ( 2C2.308), INT 

DCDF = ESTINT 

EREST = 5. 27777777 77D-6 

A ( NN ) = EREST 

GO TO 308 

IF(ERE$T.LE..5C-4)G0 TO 308 
CX = DX*(DSQRT(.25C-4/EREST) 

IF(DX.GT.DXM)GC TC 304 
DX = DXM 

X(NN4-1) = X ( NN ) 4 CX 

INT = 2 
GO TO 17C 
X( NN4-1 ) = 

IEND = 1 
GO TO 170 
CDF(NN4-1) 

NN = NN 4 



4 l.DO*. 5D0 



X ( NN ) 4CX 



4 DCDF 



910 



= CDF(NN) 

1 

IF (NN. GT .2000 ) GO TO 
GO TO ( 150.310) .IEND 
DX = DX*.0125DG/CC0F 
GO TO 160 
N = NN - 1 
NP1 = NN 
DO 400 NN=1 » N 

CDF( NN4-1 ) = CDF(NN4-1)/C0F(N4-1) 
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2 FORMAT ( ' CDF(NN) = *,F20.15) 

RETURN 

900 PRINT 1,NN»X(NN) ,X(NN+1) »EFRCR 
STOP 

1 FORMAT ( * I »XI»XI+1,ERPGR , »I5,1P3C15.6) 
910 PRINT 4 , X( NN ) 

4 FORMAT ( ' RUN CUT OF SPACE AT X=',F12.5) 
STOP 
ENO 



3 00 



350 

360 

400 

1 

410 



SUBROUTINE CDFINV(C,P) 

REAL*8 X,CDF,A, B , CC , DCC,A0INT1 , XC, E 
COMMON/CDF /X( 2001) ,CCF (2001 ) , A (2001 ), B (2001 ), NP1 
CALL TABL0C(CDF»CELE(C),NP1,I ) 

XC = X(I) + (C-CCFU ) 1/ (CDF( 1 + 1) - CDF(I))*(X( 1 + 1) - 
IF(DABS(A(I)).LE..5D-4)G0 TO 410 
CO 400 J=l,10 

OC = AOINTKX(I) ,XC,E,IEP) 

OOC = ( CDF ( I ) ♦ DC - C)/PDF(XC) 

XC = XC - DDC 

IF (DAB S(CDC) .LT..5C-4)G0 TO 410 
IF ( XC . LE .X ( I ) ) GO TO 350 
IF (XC.GE.X(I+1 ) )GC TO 360 
GO TO 4CC 

XC = X(I1 + (C - COF(I) )/DC*< XC - XC I » > 



X( ! > ) 



SP 



_ + (C - CDF(I+1) ) /(COF(I) + DC - CDF( 1+1 ) )*<X( 1+1 )-X 

CONTINUE 
URI TE ( 6 » 1 ) C , XC 
FORMAT ( • 10 ITERATIONS 
p - xf 
RETURN 
END 



IN CDF INV , C AND XC =«,2E12.4> 



SUBROUTINE TABLOC( XT, X, M,NL ) 
RE AL*8 XT ( 1 ) , X 

NT = ALGG(FLOAT(M ) /.301 + 1. 
NU = M 
NL = 1 

DO 200 1=1, NT 
NG = (NU + NL)/2 
IF ( X.GE • XT (NG) I G C TC 100 
NU = NG 
GO TO 200 
100 NL = NG 
200 CONTINUE 
RETURN 
END 



FUNCTION ADINTKXl ,X2 ,ERROR, IER) 

THIS SUBROUTINE USES DCADR5 TO OBTAIN THE INTEGRAL OF A 
PEARSON TYPE I OF II CIS T R I5UT ION FUNCTION. SUBTRACTING 
OUT THE SINGULARITY IS USED. 

IMPLICIT REAL*8 (A-Z) 

INTEGER*4 T, IER 
EXTERNAL F1,F2,PCF 

COMMON /PPARM/CO ,C1»C2»A1»A2,BC0,8C1»M1,M2»KINV»XL ,XR , MEAN ,T 
IF ( X 1 .GT .(XL+XR)/2.CQ)G0 TO 2C0 
IF ( Ml . GT .0 .DO ) GC TC 300 

PR = DCADRE(F1, XI, X2,l.C-6, 0.00, ERROR, IER) 

ADINT1 = PR + <A2-A1)**M2*( ( X2-A1 ) ** < M1+ 1. DO) - ( Xl-Al ) ** ( MH-1. DO 
1 /(Ml + l.CO)/KINV 
RETURN 

200 I F ( M2 . GT .0 .DO ) GO TC 300 

PR = DCACRE(F2,X1,X2,1. D-6 ,O.CO,ERFOR, IER) 

ADINT1 = PR- (A 2 - A1)**M1*((A2-X2)**(M2+1.00)-(A2-X1)**(M2+1.DO 
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1 / (M2+1.C0 )/KINV 
RETURN 

300 ADINT1 = 0CADRE( P0F,X1,X2, 1 . D-6, 0. CO ♦ ERROR , IEP ) 
RETURN 
END 



FUNCTION FKX) 

IMPLICIT RSAL*fi (A-Z) 

COMMON /FPARM/C0,C1,C2,A1,A2, B CO , BC1 , Ml , M2, K IN V ♦ XL , XR,MEAN,T 
IF ( X.LE «A1 )G0 TO 2 0 C 

FI = T X-Al )**M1* ( (A2 - X l **M2 - (A2 - Al> **M2) /KINV 
RETURN 

200 FI = O.DC 
RETURN 
END 



200 



FUNCTION F2( X) 

IMPLICIT REAL*8( A-Z ) 

CCMMON/P FARM/CO t Cl , C2 , A1 , A2 , 8C0, 8C 1, Ml , M2, K INV , XL , XR , MEAN , T 
IF ( X.GE «A2 IGO TO 2C0 

F2 = (A2- X )**M2*( (X-A1)**M1 - ( A 2-A 1 ) **M 1 ) /K I NV 

RETURN 

F2 = O.DC 

RETURN 

END 



200 

300 

400 



SUBROUTINE FCUTZ (FCOF,XT,NX?,FN) 

THIS SUBROUTINE GENERATES THE STATISTIC for T HS FCUTZ FN TES T . 
INPUT VARIABLES ARE : 

PCDF - THE CUMULATIVE DISTRIBUTION FUNCTION AGAINST WHICH T H? 

DEVIATES ARE BEING TESTED. CALLING SEQUENCE MUST B c OF 
THE FCRM • CALL PCCF(X,P> ', WHERE X IS AN INPU T VALUE, 
AND THE VALUE OF THE CUMULATIVE DISTRIBUTION FUNCTION 
IS RETURNED IN P. 

P MUST EE BETWEEN 0 AND 1. 

XT - THE ARRAV OF D C VIATE$, IN INCREASING ORDER. 

NXT - THE NUMBER OF DEVIATES IN THE ARRAY XT ( = N - 1) 



THE RETURNED VALUE IS FN, THE 

A 



NXT IS PRESENTLY LIMITED TO 
THE VARIABLE XTO. 

DIMENSION XT (1 ) 

REAL*8 XTD< 51) ,RN ,FND 
N = NXT + I 
00 200 I =1 ,NXT 
K * N - I 

CALL PCDF(XT(K),P) 

XTD( K+l ) = P 

RN = l.DC/N 

XTD(I) * RN - XTCC2) 

DO 300 I *2, NXT 

XTD(I) * RN - XTD ( 1+1) + XTD ( I ) 
XT D < N ) * RN - l.CO 4 XTO(N) 

FND = 0. 

CO 400 1*1, N 

FND = FND + CMAX1(XTC(I)»0.D3) 

FN = FND 

RETURN 

END 



VALUE OF THE STATISTIC. 

MAXIMUM OF 50 BY THE DIMENSION OF 
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APPENDIX 2 



CHI SQUAREC TEST 

20 PTS 30 DTS 50 PTS 





1C? 


410 


4 29 


461 


CONTROL 


5? 


2 40 


240 


229 




1? 


48 


55 


56 


PSSF 


10? 


1784 


2908 


4125 


VS 


5? 


1362 


2384 


3645 


N (0,1) 


1? 


532 


1410 


2673 







KOLMOG-SMIR 


TEST 








2 C PTS 


30 PTS 


5 C PT 




103 


491 


495 


4 63 


CCNTFOL 


53 


253 


253 


224 




13 


49 


58 


46 


PSS f 


1C3 


1636 


220 2 


3316 


V S 


5? 


1008 


1401 


2410 


N (0,1) 


13 


270 


466 


9C5 



FQUTZ F N TEST 

2 C PTS 30 PTS 50 PTS 





103 


525 


481 


488 


CCNTRCL 


5? 


225 


250 


253 




1? 


4C 


35 


62 


RSSF 


10? 


36 03 


4010 


45 9 T 


VS 


53 


2999 


3515 


4279 


N (0,1) 


1? 


1894 


2463 


346 6 



’CCNTPCL' IS TFE TEST CF INI FORM (0*1) VS UNIFCRM(C,1) 

• R S SF ' IS THE FANDOM STAEILIZED STANDARD FORM 
ALPHA = 1.00, 6ETA = C.O 



NUM8ER OF H VP C THE SI S REJECTIONS IN 5000 REPLICATIONS 
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CCMFCl 



RSS F 
VS 

N ( C ,1 ) 



CONTROL 



RSSF 

VS 

N (0,1) 



CONTROL 



RSSF 

VS 

N (0,1) 



•CONTROL* 



'RSSF* 



NIMEER OF 





CHI SQL AR ED 


TEST 






20 PTS 


30 °TS 


50 PTS 


10? 


3 SC 


421 


498 


5? 


220 


219 


242 


1? 


49 


57 


54 


1C? 


1460 


2407 


3616 


5? 


1084 


1889 


30 29 


1? 


401 


1011 


1967 





KCLMOG-SMIR 


TEST 






20 PTS 


30 PTS 


50 PTS 


10? 


488 


46 C 


500 


5? 


260 


244 


231 


1? 


64 


50 


28 


10? 


1412 


1782 


2654 


5? 


792 


1083 


1760 


1? 


222 


324 


515 





FOUTZ FA 


TEST 






20 PTS 


3 C PTS 


50 PTS 


10? 


503 


494 


5 C 1 


5? 


2 56 


257 


241 


1? 


60 


43 


42 


10? 


2863 


2272 


3962 


5? 


2145 


2678 


3422 


1? 


11 C2 


1575 


2244 



IS THE test OF UNIF0RM(0,1) VS UNIF0RN(C,1) 

IS THE RANDOM STABILIZED STANDARD FORM 
ALFFA = 1.30, EETA = 0.0 



HYPOTHESIS REJECTIONS 



IN 5000 REPLICATIONS 
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CHI SQUARFO TFS T 







20 PTS 


3C PTS 


5 C FT < 




10? 


3 95 


•ag c 


5 C 5 


CCATFCL 


53 


221 


199 


23 2 




1 % 


39 


40 


45 


RSS F 


1C? 


1441 


2116 


3266 


V S 


53 


10 54 


1598 


2692 


N ( C , 1 ) 


is 


390 


817 


1663 







KOLMOG- <MIR 


TEST 








20 PTS 


30 PTS 


50 PTS 




10? 


493 


472 


5C2 


CONTROL 


5? 


245 


215 


265 




13 


5C 


47 


52 


PS S F 


10? 


1321 


1635 


2362 


VS 


5? 


780 


991 


15C8 


N (0,1) 


1? 


210 


294 


521 







FOUTZ FA 


TEST 








20 PTS 


3 C PTS 


50 PTS 




1C? 


468 


449 


496 


CGNT RCL 


5 % 


1 88 


22 2 


237 




1? 


27 


51 


40 


RSSF 


10? 


2479 


2678 


3306 


VS 


53 


18C4 


2059 


27C5 


N (0,1) 


13 


852 


1080 


1572 



•CON TROL' 



IS THE TEST 0= UNIFQRM(0,1) VS UNIF0RM(C,1) 



* R SSF ' 



THE BANDOM STABILIZED STANDARD FORM 
ALPHA = 1.60, EETA = 0.0 



NLMEER OF HYPOTHESIS 



REJECTIONS IN 5000 REPLICATIONS 
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CHI SQUARED T ES T 
2 C °TS 30 PTS 50 PTS 





10% 


389 


4C5 


491 


CONTROL 


5 ? 


. 230 


221 


235 




1? 


43 


47 


4 € 


PSSF 


10% 


1326 


1933 


3046 


VS 


5 % 


932 


1469 


2428 


N (0,1) 


1 % 


32C 


688 


1440 







KOLMOG-S^IR 


TEST 








20 PTS 


3 C PTS 


50 PTS 




1C? 


5 27 


474 


489 


CONTROL 


5? 


275 


236 


262 




1% 


59 


38 


54 


RSSF 


10% 


1226 


1521 


2184 


vs 


5? 


760 


9 29 


1389 


t (0,1) 


1? 


242 


271 


474 



FOUTZ FN TEST 

20 PTS 30 PTS 50 PIS 





1 c* 


531 


505 


516 


CCNTFQL 


5? 


265 


263 


2 6 2 




1? 


53 


49 


54 


RSSF 


1C? 


2046 


2274 


2716 


VS 


5? 


1396 


1646 


2081 


N <C,1) 


1 % 


581 


788 


1092 



’CCNTFCL' IS 1HE TEST OF UNIF0PM(0,1) VS UNIF0PM<0,1) 

•RSSF' IS TFE RANDOM STA8ILIZEC STANDARD FOR N 
ALPHA * 1.90, SET A = 0.0 



K'PBEF OF HYPOTHESIS REJECTIONS IN 50 00 R EP L ICA TI 0N3 
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CONTROL 



RSSF 

VS 

N (0,1) 



CONTROL 



RSSF 

VS 

N (0,1) 



CONTROL 



RS SF 
VS 

N (0,1) 



'CONTROL' 



• RSSF' 



M_‘ N B E F OF 



CHI SQUARED TEST 

20 PTS 3 C PTS 50 PTS 



1C? 


4 CO 


457 


496 


5 % 


234 


242 


246 


1 % 


40 


51 


57 


10 z 


1728 


2935 


4056 


5 % 


13 16 


2435 


3629 


1% 


555 


1454 


2612 





KOLMOG— SHIR 


TEST 






20 PTS 


3 C PTS 


50 PT 


1C* 


504 


474 


477 


5 % 


252 


235 


234 


1 % 


68 


52 


47 


10* 


1565 


2222 


222 1 


5? 


966 


1410 


2372 


1* 


3 1 C 


481 


901 



FOU T Z FN TEST 
20 PTS 



ICS 


493 


53 


228 


1? 


49 


10* 


3502 


5 % 


2859 


13 


1752 



3C PTS 50 FTS 



522 


467 


255 


250 


53 


46 


4014 


4596 


2522 


42 9 9 


2547 


3487 



S THE TEST 0*= UNIFORMS, 1) VS UNIFQRM(0,1) 

IS THE RANDOM STAeiLIZEC STANDARD FCRM 
ALPHA = l.CO, BETA = 0. 25 

HYPOTHESIS REJECTIONS IN 5000 REPLICATIONS 
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CHI SCUARED TEST 



20 PTS 3 C PTS 50 PTS 





102 


372 


416 


465 


CONTROL 


52 


225 


232 


204 




12 


42 


56 


45 


RSSF 


102 


1521 


2451 


3593 


\S 


52 


1128 


1933 


3106 


N (0 ,1 ) 


12 


454 


1057 


2046 







KOLMOG-SMIR 


TEST 








20 FTS 


3 C PTS 


50 PTS 




102 


486 


512 


482 


CCNTRCL 


52 


242 


272 


2 5 E 




12 


48 


62 


45 


RSSF 


1C? 


1430 


1849 


2 7 Cl 


VS 


5? 


890 


1171 


1882 


N (0,1) 


12 


268 


271 


6 7 9 



FCUTZ FN TEST 

20 PTS 30 PTS 50 PTS 





10? 


461 


465 


455 


CONTROL 


5? 


232 


252 


226 




1? 


37 


53 


47 


RSSF 


10? 


2868 


3284 


40 26 


VS 


52 


22 36 


2684 


3486 


N (0,1) 


1? 


1174 


1545 


2325 



•CCNTRCL* IS THE T f ST OF UNI FORM ( 0 , 1) VS UNIFORM (0,1) 

•RSSF* IS TFE RANDOM STABILIZED STANDARD FORM 
ALPHA = 1.20, E ETA = C.25 



N LM E E F OF F YP C 7H E S IS REJECTIONS IN 5000 REPLICATIONS 
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CFI SQUARED 


TEST 








20 PT S 


3 C PT S 


50 PT 




1C1 


355 


390 


467 


CONTROL 


51 


223 


212 


2 42 




11 


31 


54 


53 


RSSF 


101 


14 16 


2100 


3277 


VS 


51 


1 C39 


16 10 


27C5 


N (C,l ) 


H 


357 


786 


16C7 







KOLMGG-SMIR 


T EST 








20 RTS 


30 PTS 


5 C PTS 




101 


514 


518 


5C4 


CCNTFOL 


51 


260 


282 


278 




1 1 


50 


61 


65 


RSSF 


1C1 


1345 


1661 


24C2 


VS 


5? 


807 


1 C17 


1616 


N ( C ,1 1 


11 


239 


331 


539 







FOUTZ FN 


TEST 








2 C PTS 


3 C PTS 


50 PTS 




101 


492 


5C2 


495 


CCNTPCL 


51 


248 


260 


24C 




11 


45 


35 


C - 


RSSF 


101 


23 83 


2650 


3260 


VS 


5? 


16 98 


2023 


26 1 1 


N 10,1) 


11 


822 


1027 


1473 



• CCNTRCL* IS THE TEST CF UN!FCRM(0,1) VS UNIF0RM<0,1) 

•RSSF* IS THE RANDOM STABILIZED STANDARD FORM 
ALPHA = 1. 60, BETA = C.25 



MJMeEP OF HYPCTHESIS REJECTIONS 



IN 5000 REPLICATIONS 
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CHI SQUARED TEST 







20 PTS 


30 PTS 


50 PT 




1C% 


402 


44 1 


497 


CONTROL 


5 % 


222 


230 


237 




1 % 


A3 


45 


46 


RSSF 

VS 


10 % 


1289 


1916 


3 044 


5% 


933 


1458 


2409 


N (0,1) 


11 


319 


652 


1326 



KCLNCG-SMIR TEST 



20 PT$ 3 C PTS 50 PIS 





10% 


488 


46 9 


487 


CONTROL 


CO? 


234 


235 


242 




1 % 


AO 


52 


C z 


RSSF 


10 % 


1239 


1471 


2111 


VS 


5 % 


713 


904 


1357 


N (0,1) 


1 % 


223 


269 


A A5 







FOUTZ FN 


TEST 








20 PTS 


3 C PTS 


50 PTS 




1C? 


484 


A70 


494 


CONTROL 


5? 


215 


246 


242 




1? 


56 


50 


49 


RSSF 


10% 


2003 


2289 


2773 


vs 




1358 


1688 


2098 


N (0,1) 


i% 


552 


770 


IOC? 



•CCNTRCL' IS THE TEST OF UMFCRM(O.l) VS UNIFCRM(0,1) 



•RSSF* IS THE FANDOM STAEILIZEC STANDARD FORM 
ALPHA = 1,90, EETA = 0.25 



MMEEF OF HYPOTHESIS REJECTIONS IN 5000 REPLICATIONS 
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CHI SOU AR ED TEST 

20 PTS 30 PTS 50 FT S 





10? 


3 83 


427 


455 


CC N TP CL 


5? 


213 


248 


253 




1 * 


45 


56 


37 


RSSF 


10? 


1732 


2905 


4113 


VS 


5? 


13 C3 


2395 


3668 


N (0,1) 


1? 


539 


1432 


2660 







KCLMCG-SMIR 


TEST 








20 PTS 


30 PTS 


50 PTS 




10? 


5 Cl 


49 2 


464 


CCNTFCL 


5? 


2 45 


233 


232 




1? 


38 


47 


37 


RSS F 


1C? 


1663 


2188 


3252 


VS 


5? 


1011 


1430 


2414 


N (0,1) 


1? 


289 


455 


93 C 



FOUTZ FN TEST 

2 C PTS 30 PTS 50 PTS 





10? 


507 


492 


534 


CCNTRCL 


5? 


259 


237 


271 




1? 


53 


58 


50 


RSSF 


10? 


34 52 


4011 


4571 


VS 


5? 


2399 


3512 


4297 


N (0,1) 


It 


1822 


2489 


3465 



•CCNTRCl' IS THE TEST CF (JNI FORM <0 , 1 ) VS UNIFCP.M(C,I) 

•RSSF' IS THE FANDOM STAEILIZED STANDARD FORM 
ALPHA = I. CO, 0ETA = 0.50 



NUNS ER OF HYPOTHECS REJECTIONS IN 5000 REPLICATIONS 
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CHI SOLAR ED TEST 



20 PTS 3 C PTS 50 RTS 





1 c - 


382 


406 


443 


CONTROL 


5 % 


213 


19 5 


2 C 7 




1 % 


45 


38 


3 6 


RSSF 


10 % 


1665 


2552 


36S5 


V S 


5 % 


1266 


2038 


3222 


N (0,1) 


1 % 


509 


i 178 


2.2 02 







KCLMCG-SMIR 


TEST 








20 PTS 


3 0 PTS 


50 PTS 




10? 


501 


5 03 


486 


CONTROL 


5? 


248 


256 


2 53 




1? 


45 


55 


57 


PSSF 


10? 


16CC 


2092 


30 2 5 


VS 


5? 


1026 


1425 


2208 


N (0,1) 


1? 


342 


532 


951 



FOUTZ FN TEST 

20 PT S 3 C PTS 50 °TS 





1C? 


5 06 


467 


471 


CONTROL 


5? 


246 


236 


228 




1? 


44 


46 


43 


RSSF 


10? 


28 16 


2278 


3920 


VS 


5 % 


2168 


2680 


3382 


N (0,1 ) 


1 % 


1180 


1546 


2270 



•CCMPCL* IS 11- E TEST OF UMFCRM<0,1) VS UNIFORMS, 1) 

' RSSF* IS THE PANDCM STAeiLIZEC STANDARD FORM 
ALPHA = 1.30, EETA = 0. 50 



M.PBEF CF HYPOTHESIS REJECTIONS IN 5000 REPLICATIONS 
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CHI SQUARED TEST 



CONTROL 



RSSF 

VS 

N (0,1) 



CONTROL 



PSSF 

VS 

N (0,1) 



CONTROL 



RSSf 

VS 

N (0,1) 



•CONTROL' 



• RSS F* 



NLMEER OF 





20 PTS 


3 C PTS 


50 PTS 


1C? 


413 


440 


470 


5 ? 


251 


222 


2C5 


1 % 


36 


48 


41 


1 C? 


1451 


2269 


3435 


5? 


1063 


1794 


2841 


1? 


366 


948 


1626 





KCLNCG-SMI R 


TEST 






20 PTS 


3 C PTS 


50 PTS 


10? 


478 


518 


513 


5? 


2 52 


278 


277 


1? 


49 


65 


6C 


10? 


1435 


1853 


2636 


5? 


888 


1229 


1881 


1? 


2 67 


444 


758 





FOUTZ FN 


TEST 






20 PTS 


3 C PTS 


50 PTS 


1C? 


493 


463 


488 


5? 


241 


220 


227 


1? 


47 


47 


40 


10? 


2314 


2643 


3268 


5? 


16 86 


202 8 


2594 


1? 


761 


1 C52 


1449 



S THE TEST OF UMFCRM(O.l) VS J NI FORM (0, 1 ) 



IS THE RANDOM ST A El LI 
ALPHA = 1.60, EETA 



HYPOTFES IS REJECTIONS 



EC STANDARD FORM 
= 0*50 



IN 5000 R EP LI CAT I C NS 
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CHI 5QU AR EO TEST 



20 PTS 30 PTS 5 C PTS 





1C? 


419 


411 


466 


CCNTFOL 


5? 


218 


230 


234 




135 


38 


54 


52 


RSSF 


1C? 


1299 


1934 


3 C42 


VS 


5? 


9 16 


1438 


24C7 


N ( C,1 ) 


1 % 


297 


680 


1348 







KOLMOG- SMIR 


TEST 








20 PTS 


30 PTS 


50 PTS 




10% 


5 C2 


494 


48 2 


CONTPCL 


5% 


26C 


255 


248 




1? 


47 


58 


62 


RSSf 


1C? 


1239 


1559 


2146 


VS 


5? 


738 


961 


1357 


N (C,l) 


1? 


214 


296 


475 







FCUTZ FN 


TEST 








2 0 PTS 


30 FTS 


50 PTS 




10% 


473 


479 


508 


CONTROL 


5% 


259 


229 


240 




1% 


5 5 


53 


52 


RSSF 


10? 


2018 


2310 


2708 


VS 


5? 


1384 


1629 


2043 


N (0,1) 


1? 


591 


749 


1023 



•CONTROL' IS THE TEST OF UNIFORM(0,1> VS UNIFORMOtl) 

* R S SF* IS THE FANDOM STAEILIZEO STANDARD FORM 
ALPHA = i.90, E ETA = 0.50 



MM8EF OF FYPCTHESIS REJECTIONS IN 5000 REPLICATIONS 
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CHI SQUARED TEST 



20 PTS 3 C PTS 50 PTS 





1C? 


A 10 


445 


436 


CONTROL 


c ? 


229 


236 


221 




1? 


30 


49 


46 


RSSF 


10? 


1831 


2913 


4 1 C 8 


VS 


5 % 


13 58 


239 1 


3655 


b (C,li 


1 % 


566 


* 93 8 


2690 







KOLMOG-SMIR 


T EST 








2 C FT S 


30 PTS 


5 C PTS 




10? 


524 


505 


527 


CCNTfQL 


5? 


250 


259 


25 I 




1? 


46 


53 


55 


PSSf 


1C? 


1627 


2196 


33C4 


V s 


5 % 


10 51 


1421 


2387 


N 10,1) 


1? 


3 04 


500 


931 



F CUT Z FN TEST 



2 C PTS 30 PTS 50 PTS 





10? 


515 


497 


502 


CONTPCL 


5? 


229 


254 


2 c c 




1 % 


51 


59 


46 


PSSF 


10 % 


3586 


4057 


4542 


VS 


5 % 


3005 


2591 


4256 


N <0,i) 


1 % 


1893 


2495 


3453 



* CCNTRCL' IS THE T£$T CF UNIFORMS, 1) VS UNlFORM(Otl) 



' R S SF * 



IS THE RANDOM STABILIZED STANDARD FORM 
ALPH* = 1.00, EETA = C.75 



NUMBER OF 



HTPCTHESIS REJECTIONS 



IN 5000 REPLICATIONS 
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CHI SOUAREO TEST 



20 PTS 2C PTS 50 PTS 





1C? 


368 


447 


505 


CCNTPOL 


5 % 


2 CC 


241 


222 




1? 


43 


58 


39 


PSSF 


101 


19 41 


2511 


4 C97 


VS 


51 


1495 


2429 


36 T C 


N ( C,l) 


1? 


704 


1535 


2762 







KQLMOG-SMIR 


TEST 








20 PTS 


3C PTS 


50 PTS 




10? 


4 75 


50 4 


524 


CCNTPCL 


5? 


242 


237 


264 




1? 


47 


45 


47 


RSSP 


10? 


20 17 


2 64C 


3726 


VS 


51 


1331 


1930 


3089 


N <C,1) 


1? 


546 


904 


1762 



FCUTZ FN TEST 

20 PTS 30 PTS 50 PTS 





10? 


461 


513 


490 


CCNTPOL 


5 1 


226 


248 


248 




1 % 


45 


54 


5 5 


RSSF 


1C? 


2895 


2212 


4 C 5 1 


VS 


5? 


2227 


2712 


3542 


N (0 ,1 ) 


11 


12C9 


1595 


2426 



'CCNTPCL * IS THE TEST CF UNIFORM(Otl) VS UNIFCPMC,1) 

•PSSF* IS THE RANDOM STAEILIZEC STANDARD FORM 
ALPHA = 1.20, EET A = C.75 



HUMBER CF HYPOTHESIS REJECTIONS IN 5GC0 REPLICATIONS 
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CHI SQUARED TEST 

20 PTS 3 C PTS 50 PTS 





1C? 


4 CO 


407 


501 


CONTROL 


5? 


2 43 


225 


256 




IS 


48 


51 


46 


RS SF 


10? 


1628 


2451 


3633 


VS 


5? 


1240 


15 6 5 


3C99 


N (0,1) 


1 ? 


509 


1126 


21 1 C 







KOLMOG-SNIR 


TOST 








20 PTS 


3 C PTS 


50 FTS 




1C? 


491 


505 


476 


CCNT FOL 




111 


244 


243 




1? 


48 


54 


50 


F S SF 


10? 


1622 


2 IC4 


302 2 


VS 


5? 


10 50 


1465 


2330 


b (0,1) 


1? 


375 


588 


1145 



FOUTZ FN TEST 

20 PTS 30 PTS 50 FTS 





1C? 


522 


490 


45 C 


CONTROL 


5? 


266 


26 7 


243 




1? 


52 


60 


56 


RS SF 


10? 


2426 


2722 


3367 


VS 


5? 


1753 


20 67 


2727 


N (0,1) 


1? 


8 25 


1C58 


1614 



•CONTROL* IS THE TEST OF UNIFORM(Otl) VS UNIFORM(Otl) 

•RSSF* IS THE RANDOM STAEILIZEC STANDARD FORM 
ALPHA = 1 . 60 1 3ETA = 0.75 



M'NBEF OF HYPOTHESIS REJECTIONS IN 5000 R?d L I CAT I C HS 



CHI SCU6RED TEST 









2 C PTS 




30 PTS 


50 PT 




10* 




405 




4C7 


469 


CONTROL 


5* 




?38 




229 


215 




1* 




54 




49 


46 


FSSF 


10* 




1366 




1989 


3088 


VS 


5* 




987 




1450 


2487 


N 10,1) 


1* 




341 




704 


141 S 








KOLVOG- 


SMIR 


TEST 










20 PTS 




3 C PTS 


50 PT 


CONTROL 


10* 




50 1 




505 


483 


5* 




246 




268 


245 




1* 




57 




53 


46 


PSSF 


10* 




1273 




1558 


2246 


VS 


5* 




797 




936 


1465 


N (0,1) 


1* 




234 




308 


508 








FOUTZ F 


N TEST 










20 PTS 




3C PTS 


50 FT 




1C* 




497 




457 


478 


CONTROL 


5* 




242 




227 


242 




1* 




58 




48 


48 


PSSF 


10* 




2102 




2256 


2745 


VS 


5* 




1479 




16 15 


2 C E 1 


N (C,l) 


1* 




647 




76 0 


1076 

i 


;CNTROL* IS THE 


TEST 


OF UNI 


FCRM< 


:o,i) VS U 


NI FORM (0 


• RSS F' 


IS THE 


RANDOM STABILIZED STANDARD FCRM 




ALPHA = 


1 » SC, 


EETA 


= 0.75 




'LM9EF OF 


HYPCTH 


ESIS 


REJECT 


IONS 


IN 5000 R 


EPLICATI 



63 



CHI SQUARED TEST 



20 PTS 3 C PTS 50 PTS 





1C? 


3 55 


441 


455 


CONTROL 


5% 


213 


238 


2 43 




IS 


35 


53 


62 


PSSF 


10? 


1793 


2 52 C 


4109 


VS 


5? 


1355 


2 38 6 


3658 


b (0*11 


1? 


554 


1413 


2707 







KGLMOG-SMIR 


TEST 








20 PTS 


20 P^S 


50 FT 




10? 


493 


504 


457 


CCNTFCL 


5% 


241 


231 


228 


1% 


50 


42 


43 


RSSF 


ICS 


1640 


2203 


3325 


VS 


5% 


10 04 


139 2 


2434 


N ( C ,1 ) 


IS 


295 


421 


935 







FOUTZ FN 


TEST 








20 PTS 


30 PTS 


50 PT 




10? 


526 


488 


48 1 


CCNTRCL 


5? 


260 


247 


223 




I? 


54 


58 


50 


RSS f 


1C? 


3530 


4001 


4577 


VS 


5? 


29 33 


3501 


4284 


N ( C »1) 


IS 


1849 


2405 


34C5 



•CCNTRCL' IS THE TEST OF UNIFORM(0,1) VS UNI FORM (0 , i ) 

• FSSF • IS THE RANDOM STABILIZED STANCARD FORM 
ALPHA = 1 .00 t BETA = 1.00 



NLME EP QF HVPCTHESIS REJECTIONS IN 5000 REPLICATIONS 
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CHI S QU AR £0 TEST 

20 P T S 2C PTS 50 PTS 





1 Cl 


363 


428 


518 


CONTROL 


5? 


218 


226 


242 




II 


45 


44 


55 


RSSF 


1C? 


2651 


2826 


4682 


VS 


51 


2291 


3457 


4472 


N ( C ,1 ) 


1? 


1344 


2572 


3935 







KOLMOG-SMIR 


T EST 








20 PTS 


30 PTS 


50 PTS 




10? 


526 


533 


5 24 


CONTROL 


5? 


2 43 


263 


2 62 




1? 


51 


48 


1C 


PSS f 


10? 


2937 


2776 


4636 


vs 


5? 


23 26 


3235 


4378 


N <C,1) 


1? 


1245 


2035 


3452 



F CUT Z F N TEST 

20 PTS 30 PTS 50 PTS 





10? 


488 


477 


500 


CONTROL 


5? 


259 


216 


236 




1? 


54 


53 


44 


PSS F 


10? 


2962 


3429 


4162 


VS 


5? 


2260 


2813 


3660 


N tO,l) 


1? 


1237 


170 5 


25 50 



'CONTROL * IS THE TEST OF UNIF0RM(0,1) VS UNIFORMS, 1) 

* RSSF' IS TFE FANDOM STAEILIZEC STANDARD FORM 
ALPHA = 1. 20, E ETA = 1 *00 



MMEEP OF HYPOTHESIS REJECTIONS IN 5000 REPLICATIONS 
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CCNT FOL 



RSS f 
VS 

N < C ,1 ) 



CCNT PG L 



RSSF 

VS 

N (0,1) 



CCNTPCL 



RSSF 

VS 

N (C,l) 



• CONTROL ' 



•RSSF 



NCMEER 0 



CHI SQU AR ED TEST 

20 PTS 3C PTS 50 FT $ 



101 


A 12 


443 


472 


51 


220 


237 


220 


11 


30 


50 


46 


1 Cl 


1868 


2768 


3 S 5 8 


51 


1487 


2294 


3541 


11 


701 


1427 


2578 





KCLMOG-S 


MIR TEST 






20 PTS 


30 PTS 


50 PT 


101 


513 


477 


5 14 


51 


268 


251 


251 


11 


56 


56 


49 


101 


1968 


2603 


3652 


51 


1370 


1903 


30 16 


11 


586 


905 


1765 




FOUTZ F fi 


TEST 






20 PTS 


3 C PTS 


50 PT 


101 


5 07 


504 


505 


51 


250 


240 


252 


11 


58 


45 


58 


101 


2414 


2792 


3409 


51 


1786 


2120 


2771 


11 


853 


1123 


1625 



IS THE TEST CF UMIFORM(Otl) VS UNIF0RM(C,1) 

IS THE RAN COM STABILIZED STANDARD FO«M 
ALPHA = 1.60, EETA = 1 .00 



HYPOTHESIS REJECTIONS 



IN 5000 REPLICATIONS 
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CHI SQLARED 


TEST 








20 FTS 


3 C PTS 


50 PT 




ic* 


3 74 


431 


525 


CONTROL 


5? 


231 


232 


2 6 3 




1 % 


48 


57 


58 


FSSF 


10 % 


1329 


15 53 


3 1C 7 


VS 


5 % 


977 


1438 


2466 


N (0,1) 


1 % 


347 


708 


1459 







KQLMOG-SMIR 


TEST 








20 PTS 


3 C PTS 


50 FTS 




1C? 


494 


48 1 


516 


CCNTPOL 


5 % 


?-3g 


240 


241 




1 % 


"54 


47 


43 


RSSF 


10? 


1322 


1625 


23 C c 


VS 


5? 


8C8 


978 


1533 


N (C,U 


1 % 


225 


3 02 


540 







FOUTZ FN 


T EST 








20 FTS 


30 PTS 


5 C PTS 




10? 


525 


523 


520 


CCNTFCL 


5? 


236 


271 


296 




1? 


50 


67 


55 


RSSF 


1C? 


2086 


2 280 


2758 


VS 


5 % 


1409 


1645 


2075 


N ( C ,1 ) 


1? 


620 


71 5 


1078 



*C CNT F CL ' IS THE TEST OF UNIFQRM(C,1) VS UMPCPM(0,1) 

'FSSF ' IS THE FANCCM STABILIZEC STANDARD FORH 
ALPHA = 1 .90 t EETA = 1 .00 



N L H 6 E P CF HYPOTHESIS REJECTIONS IN 50CO REPLICATIONS 
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CHI SQUARED TEST 



CCNTPGL 
N IXNCRV 

vs 

N ( C ,1 ) 



CONTROL 



MIXNCRM 

VS 

N (C,l) 



CCNTROL 



MXNCPM 

VS 

N (0,11 





2 C PTS 


3 G PTS 


50 PTS 


10? 


388 


469 


4 £6 


5? 


2 29 


238 


224 


1? 


39 


48 


43 


1C? 


1272 


1925 


30C4 


5? 


9 09 


1407 


2381 


1? 


3 29 


664 


1337 





KOLMOG-SMIR 


TEST 






20 PTS 


30 PTS 


50 PTS 


10? 


52 5 


513 


5 10 


5? 


2 54 


260 


2 4 S 


1? 


42 


46 


64 


10? 


1235 


1537 


2133 


5? 


7 13 


S3 1 


1369 


1? 


208 


28 5 


44 S 





FOUTZ FN 


TEST 






20 PTS 


3 C PTS 


50 PTS 


10? 


508 


476 


523 


5? 


245 


222 


255 


1? 


50 


38 


6C 


10? 


1982 


2191 


26 10 


5? 


1336 


1576 


1955 


1? 


533 


672 


977 



•CCNTROL* IS THE TEST CF UNIFQRM(0,1) VS UNIFORMOtl) 

•MXNCRM IS TEE MIXED NORMAL DISTRIBUTION 
C. C *N (.0, 1) + 1 .OC*N< C,2 ) 



rUNBEP OF HYPOTHESIS REJECTIONS IN 5000 REPLICATION 



68 



CHI SQUARED TEST 



20 PTS 30 PTS 5 C FTS 





10? 


4CS 


All 


458 


CCNTFCL 


5? 


233 


205 


234 




1? 


29 


34 


At 


f'lXNCRM 


10? 


2513 


3558 


46 2 C 


VS 


5* 


2088 


3C8 0 


4342 


N (0,1) 


12 


1028 


2052 


36 C 5 







KOL NOG- SMI R 


TEST 








20 PTS 


3 C PTS 


50 PTS 




10? 


475 


463 


50C 


CONTROL 


5? 


238 


229 


242 




1? 


46 


35 


38 


NIXNCPM 


10? 


21 88 


2827 


39 6 5 


VS 


5 ? 


1420 


1974 


3161 


N (0*1) 


1? 


512 


727 


1515 



FOUTZ F 1^ TEST 

20 PTS 3 C PTS 50 PTS 





10? 


522 


454 


481 


CONTROL 


5? 


236 


217 


226 




1? 


A7 


41 


52 


MXNCFM 


10? 


3338 


3712 


4346 


VS 


5? 


2724 


3174 


3935 


N (0,1) 


1? 


1657 


2064 


2970 



' CONTROL' IS THE TEST CF UNIFORMS, 1) VS UN!F0RM(C,1) 

•MXNCRM* IS THE MIXEC NORMAL DISTRIBUTICN 
C.C *M.0,1) ♦ 1.00*N<0,3) 



MJNBEF OF HYPOTHESIS REJECTIONS IN 5000 R EP L IC A TIONS 
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CHI SQUARED TEST 

20PTS 3CPTS 50PTS 





10? 


407 


425 


5 11 


CONTROL 


5 % 


2 34 


232 


252 


1 % 


48 


4 £ 


5C 


MIXNOFM 


10? 


3427 


4423 


4925 


VS 


5 % 


3006 


4154 


4863 


N (0 ,1) 


1? 


1 7 £1 


3377 


46 29 







KOLMOG-SMIR 


TEST 








20 PTS 


3C PTS 


50 FTS 




1C? 


510 


490 


511 


CCNTFOL 


5 % 


261 


241 


250 




1? 


59 


48 


47 


MXNCRM 


1C? 


2976 


2799 


4717 


VS 


5? 


2113 


3025 


42 £6 


N ( Ctl ) 


1? 


890 


1471 


2792 



FOUTZ FN TEST 

20 FT S 30 PTS 5C PTS 





10? 


518 


501 


528 


CONTROL 


5 % 


242 


268 


276 




1 % 


68 


66 


62 


MIXNCRM 


1C? 


4131 


4502 


4 £6 5 


VS 


5? 


2686 


4195 


4742 


N (C,l) 


1? 


2653 


2330 


4272 



•CCNTRCC* IS THE TEST OF UN I^ORM ( C ,1 ) VS UNIFCRM<C,1> 

•MIXNCRM' IS THE MIXED NORMAL DISTRIBUTION 
0.0 *N< .0»1 ) + 1 »00*N ( 0 »4 1 



NUMBER OF HTPCTHESIS REJECTIONS IN 5000 REPLICATIONS 
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CHI SOU £R ED 


TEST 








20 PTS 


20 PTS 


5 C PTS 




10? 


3 82 


422 


510 


CCNT POL 


5? 


215 


137 


254 




IZ 


41 


38 


46 


HmCRH 


10? 


472 


542 


713 


vs 


5% 


2 53 


308 


3 70 


F (Ctl) 


1% 


29 


66 


74 







KOLMOG-SMI R 


TEST 








20 PTS 


30 PTS 


50 PTS 




1C? 


488 


537 


525 


CCNTFCL 


5% 


238 


250 


265 




1? 


46 


50 


53 


M IXNCRM 


1C? 


587 


667 


655 


VS 


5? 


287 


328 


367 


h ( c a ) 


1? 


54 


70 


84 







FOUTZ FN 


TEST 








20 FT S 


30 PTS 


50 PTS 




10? 


468 


480 


517 


CCIVTPCI 


5? 


220 


258 


247 




1? 


46 


64 


56 


MIXNCRN 


1 C? 


799 


796 


8 El 


VS 


5? 


425 


426 


477 


N (oai 


1? 


105 


107 


136 



• C CNTPC L ' IS THE TEST CF UNIFORMS, 1) VS UNIF0RM(C,1) 

' M IXNQRM * IS THE MIXED IVCPMAL DISTRIEUTION 
C.7C*N(.0,1> + C.30*N(0,2> 



M.MEEP OF HYPOTHESIS REJECTIONS IN 5000 REPLICATICNS 
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CHI SQUARED TEST 



20 PTS 3C PTS 5 C PTS 





10 ? 


421 


457 


510 


CCNT FOL 


5? 


245 


244 


22 1 




1% 


49 


52 


3 € 


MIXNCRM 


1C? 


551 


727 


1C4 7 


VS 


5? 


345 


438 


602 


N ( C ,1 ) 


1 % 


82 


131 


185 







KOLMCG-SMIR 


TEST 








20 PTS 


30 PTS 


50 PTS 




10? 


490 


504 


527 


CCNTROL 


5? 


251 


256 


23 S 




1? 


46 


58 


42 


MIXNORM 


10? 


6 46 


745 


845 


VS 


5? 


337 


405 


447 


N (0.1) 


1? 


81 


93 


E3 



FQUTZ F I' TEST 

20 PTS 3 C PTS 50 PTS 





10? 


5 19 


494 


507 


CONTROL 


5? 


24C 


242 


256 




1? 


55 


46 


51 


MIXNCRM 


10? 


1029 


1139 


1249 


VS 


5? 


59 1 


671 


778 


h (0,1) 


1? 


172 


193 


262 



•CCNTRCL* IS TFE TEST CF UNIFORMS, 1) VS UNI FOR V ( Ct I ) 

•MIXNORM IS THE MIXED NORMAL DISTRIBUTION 
0.7 C*N (.0,1) + 0.30*N(0»3) 

KIMBEF CF HYPCTHESIS REJECTIONS IN 5000 REPLICATIONS 
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CHI SOU JR EO TEST 



20 PTS 3 C P7S 50 PTS 





1C? 


400 


435 


5 15 


CONTROL 


5? 


237 


223 


240 




1? 


45 


50 


52 


M I > N CR M 


10? 


6 36 


910 


1320 


VS 


1 % 


4C5 


564 


336 


N (0,1 ) 


1 % 


109 


175 


222 







KOLMOG-SMIR 


TEST 








2 C PTS 


30 P T $ 


50 PT 




10? 


508 


531 


5 C E 


CCNTFCL 


5 % 


259 


259 


248 




1 % 


54 


51 


54 


N IXNCRM 


1C? 


747 


853 


10C1 


VS 


5? 


407 


441 


532 


N (C,l ) 


1? 


1 15 


122 


148 



FO'JTZ FN TEST 

20 FT 5 30 PTS 50 PTS 





10? 


4 92 


478 


511 


CONTROL 


5 % 


2 29 


232 


241 




1? 


61 


47 


42 


MIXNCRM 


10? 


1281 


1418 


166E 


vs 


5? 


774 


888 


1121 


N (0,1) 


1? 


238 


222 


4 C 1 



* CCNTRC L* IS THE TEST OF UNI FOR M ( 0 ♦ 1 ) VS UNIFORMS, 1) 

' M I XNORM • IS THE MIXED NCRMAL DISTRIBUTION 
C.7C*N(.0,1) + 0.30*N<0,4) 



MMEEP OF HYPOTHESIS REJECTIONS IN 5000 REPLICATIONS 
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CHI SQU £R ED TEST 







20 PTS 


30 PTS 


50 FT S 




10? 


3 80 


431 


5C7 


CCNTFCL 


53 


221 


245 


220 


13 


40 


59 


56 


MXNCRM 


1C? 


388 


484 


545 


V S 


5? 


214 


274 


275 


N <G,1) 


1? 


40 


63 


62 






KOLMOG- 


- SM IR TEST 








20 PTS 


30 PTS 


50 PTS 




10? 


468 


480 


48 1 


CCNTRCL 


5? 


220 


236 


232 




1? 


53 


46 


28 


MIXNCRM 


10? 


519 


539 


556 


VS 


5? 


258 


277 


286 


N (0,1) 


1? 


58 


55 


53 






FOUTZ 


Ft TEST 








20 PTS 


3 C ° T S 


50 PTS 




10? 


487 


498 


502 


CONTROL 


5? 


234 


242 


265 




1? 


42 


45 


54 


VIXNC FM 


10? 


684 


689 


687 


VS 


5% 


352 


354 


400 


N (0,1) 


1? 


83 


92 


96 



•CONTROL* IS THE TEST CF UNIP0RM(0,1) VS UNIFORM(C,l) 

•NIXNCPM* IS THE MIXED NORMAL DISTRIBUTION 
C. €C*N ( .0* I ) + 0.20*N(0,2) 



t'UMBEF OF HYPOTHESIS REJECTIONS IN 5000 REPLICATICNS 
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CHI SQUARED TEST 

20 PTS 30 PTS 5 C PTS 





1C? 


417 


405 


*92 


CONTROL 


5? 


218 


211 


242 




1? 


43 


38 


42 


b IXNCRM 


1C? 


477 


544 


736 


VS 


5? 


281 


306 


364 


N ( C,1 ) 


1% 


48 


64 


58 







KOLMOG-SMIR 


T EST 








20 FTS 


30 P T S 


50 PT 


CCNTRCL 


10? 


5 CC 


469 


527 


5? 


267 


222 


2 6 C 




1? 


47 


46 


48 


MIXNCRM 


1C? 


587 


590 


678 


vs 


5? 


3 17 


297 


342 


N * ( C tl ) 


1? 


69 


59 


78 



FCUTZ FN TEST 

20 PTS 30 PTS 50 PTS 





10? 


495 


46? 


479 


CONTROL 


5? 


234 


213 


245 




1? 


49 


44 


48 


NIX NC PM 


10? 


839 


819 


945 


VS 


5? 


475 


454 


522 


N (Ctl) 


1? 


1 18 


130 


142 



’CCNTRCL ' IS THE TEST OF UNIF0RM(0»1) VS UN IFQR N < 0 , 1 ) 

•MIXNCRN' IS THE MIXED NORMAL DISTRIBUTION 
C • £C*N <• Ot 1) + 0.2C*N<0,2) 



NUMBER OF HYPOTHESIS REJECTIONS IN 5000 REPLICATIONS 
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CHI SQUARED TEST 







20 PTS 


3 C PTS 


50 PTS 




1C* 


402 


4 50 


539 


CCMFCL 


5? 


225 


242 


247 


1* 


43 


49 


6C 


MIXNGRM 


1C* 


5 07 


637 


9C2 


\S 


51 


310 


361 


496 


N <C,1) 


11 


74 


89 


138 






KOLMOG- SMIR 


TEST 








20 PTS 


30 PTS 


50 PTS 




10? 


511 


502 


542 


CONTROL 


51 


253 


252 


270 




1? 


51 


43 


63 


NIXNCRM 


10? 


622 


661 


712 


VS 


5? 


3 29 


360 


416 


N (0,1) 


1? 


85 


69 


93 







FCUTZ FN 


TEST 








20 PTS 


30 PTS 


50 PTS 




10? 


485 


496 


514 


CONTROL 


5? 


238 


243 


259 




1? 


48 


48 


c c 
<0* +* 


M IXNQRM 


10? 


958 


1011 


1168 


VS 


5 % 


548 


583 


721 


N (0,1) 


1? 


171 


174 


246 



•CONTROL* IS THE 



TEST CF UNI FCRM (0*1) VS UNIFORMO,!) 



• M 1 >NCR V’ IS TFE MIXED 
C.£C*Nl.O,l) + 



NORMAL D 1ST RI8UTI ON 
0.2 0*N< C,4 ) 



H NB E F CF HYPOTHESIS REJECTIONS IN 5000 REPLICATION 
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CHI SQUARED TEST 

20 PTS 3 C FT S 50 PTS 





1C3 


420 


470 


47C 


control 


53 


231 


240 


236 




1? 


40 


51 


44 


MI >NCPM 


10? 


4 1 C 


478 


525 


VS 


53 


228 


241 


247 


N (0,1) 


13 


35 


48 


46 







KOLMOG-SNIR 


TEST 








20 PTS 


3 C PTS 


50 PTS 


CONTROL 


1C? 


499 


490 


519 


53 


244 


245 


251 




13 


54 


44 


55 


MXNCRM 


10? 


521 


512 


535 


VS 


5? 


264 


257 


266 


N <C,1> 


1 3 


57 


47 


64 



FOUTZ FN TEST 

20 PTS 30 PTS 50 PTS 





1C? 


493 


490 


4 £ C 


CONTROL 


53 


245 


256 


229 




13 


43 


62 


47 


HXNCRN 


10? 


577 


, 583 


587 


\S 


53 


282 


310 


2 S C 


N (C,l) 


13 


66 


73 


77 



•CONTROL* IS THE TEST OF UNIFCRMO.l) VS JNIFQRMOtl) 

•MIXNORM* IS THE MIXED NORMAL E ISTP. IBUT I CN 
0.907NI .0,1.) + 0.10*N(0,2) 



NUME EP OF HYPOTHESIS REJECTIONS IN 5000 REPLICATIONS 
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CHI SQU £R ED TEST 

20 PTS 30 PTS 50 FT S 





10? 


432 


416 


4 £6 


CONTROL 


5? 


246 


214 


215 


1? 


43 


41 


4 C 


V IXNCRM 


10? 


4 56 


431 


537 


VS 


53! 


262 


233 


2 58 


n cc,n 


1? 


51 


43 


53 







KOLMOG-SMI R 


T EST 








20 PTS 


30 PTS 


50 PTS 




10* 


5 C9 


495 


509 


CCNTFOL 


5? 


271 


239 


25 2 




1? 


55 


46 


55 


MIXNCRM 


1C? 


553 


5-26 


557 


VS 


5? 


302 


26 1 


274 


N (C ,1 ) 


1? 


67 


53 


64 



FOUTZ FN TEST 

20 PTS 3 C PTS 50 PTS 





10? 


516 


471 


520 


CCNTFOL 


5? 


271 


239 


252 




1? 


58 


51 


47 


MIXNCRM 


10? 


66 C 


653 


675 


VS 


5? 


353 


357 


356 


N (0,1) 


1? 


95 


88 


£ 5 



•CONTROL' IS THE TEST CF UNIF0RM<0, 1) VS UNIFCRM(0,1J 

•NI>NCPy IS THE MIX EC NORMAL DISTRIBUTION 
C*5C*N(.0,1) + 0.1C*N(0,3) 



NUMEER OF HYPOTHESIS REJECTIONS IN 5000 REPLICATIONS 
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CHI SQUARED TEST 

20 PTS 30 PTS 50 PTS 





1C? 


357 


439 


505 


CON TROL 


5? 


240 


244 


246 




1? 


31 


53 


55 


R IXNCRN 


1C? 


419 


50 8 


551 


VS 


5? 


260 


266 


3C2 


N (C,l) 


1? 


36 


65 


74 







KCLMQG-SMIR 


T EST 








2 C PTS 


3C PTS 


5 C PT 


CCNTfOL 


10? 


4 84 


516 


48 £ 


5? 


240 


247 


229 




1? 


56 


47 


45 


HIXNCRN 


1C? 


525 


534 


54 £ 


V s 


5? 


268 


28 2 


27 7 


N (0,11 


1? 


64 


61 


6 C 



FO'JTZ FN TEST 



2 C PTS 30 PTS 50 PTS 





10? 


514 


511 


45 1 


CCNTRCL 


5? 


7 47 


253 


245 




1? 


48 


50 


58 


NIXNCRM 


10? 


736 


707 


747 


VS 


5? 


385 


384 


422 


N (0,1) 


1? 


81 


105 


112 



•CCNTRCL' IS THE TEST CF UNIFORMS, 1) VS UNIFCRM(C,1) 

•RIXNCRN' IS THE HIXEC NCRMAL CISTRIBUTION 
C. 90*N ( .0,1 ) + 0.10*N<0,4) 



NUMBER CF HYPOTHESIS REJECTIONS IN 5000 REPLICATIONS 
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CHI SQUARED 


TEST 








20 PTS 


30 PTS 


50 PTS 




10% 


3 €6 


417 


494 


control 


5% 


225 


241 


238 




1% 


39 


44 


46 


R IXNCRR 


10% 


1459 


2008 


3146 


VS 


5% 


1 093 


15 3 3 


2556 


N (C,l) 


1% 


407 


793 


1568 







KCLMOG-SMIR 


TEST 








20 PTS 


3C PTS 


5 C PTS 




10% 


493 


483 


487 


CCNTFOL 


5% 


267 


246 


238 




1% 


52 


53 


36 


MIXNCRR 


1C% 


2179 


2806 


3756 


VS 


5% 


16C6 


2181 


3226 


N <C,1) 


1% 


714 


1148 


2 043 







FOUTZ FN 


TEST 








2 C PTS 


30 PTS 


50 PTS 




10% 


535 


495 


499 


CONTROL 


5% 


250 


246 


248 




1% 


54 


43 


41 


MIXNGRM 


10% 


18 09 


1994 


2492 


VS 


5% 


1198 


1428 


1825 


N (0,1) 


1% 


509 


60 2 


865 



’CCNTRCL' IS THE TEST CF UNIFORMED ,11 VS UNIFORM (0,1) 

* M I >NO PM ' IS THE MIXEO NORMAL DISTRIBUTION 
C» 7 C*N ( .5, 1) + 0 .3C*N( C,3) 



NUMBER OF HYPOTHESIS REJECTIONS IN 5000 REP LICATI CHS 
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CHI SQUARED TEST 

20 PTS 30 PTS 50 FT S 





102 


377 


449 


512 


CCNTPCL 


52 


20 5 


239 


257 




12 


46 


46 


44 


NIXNCRM 


1C? 


1557 


2112 


33C4 


VS 


5? 


1170 


1563 


2660 


N 10,1) 


1? 


473 


750 


1 6 4 C 







KOLf'CG-SMI R 


TEST 








20 PTS 


30 PTS 


50 PTS 




10? 


482 


481 


538 


CONTROL 


5? 


229 


260 


271 




12 


45 


69 


56 


MIXNCRM 


10? 


24 5 C 


3101 


4056 


VS 


5? 


1855 


2510 


3596 


N (0,1) 


1? 


9 Cl 


1344 


24CQ 



FOUTZ FN TEST 

20 PTS 3 C PTS 50 PTS 





1C? 


479 


512 


531 


CONTROL 


5? 


233 


262 


255 




12 


50 


53 


53 


MXNCRN 


10? 


1781 


1917 


2405 


VS 


5? 


12C9 


1369 


1774 


N (0,1) 


1? 


462 


573 


842 



•CCNTRCL' IS THE TEST CF UNI FORM (0 » 1) VS UN I FCR M ( C , I ) 

• ^ I X N CR N • IS THE MIXEC NORMAL CISTP.IBUTICN 
0.£C*M.5,1> + 0 . 20* N( 0,3) 



NUMBER CF HYPOTHESIS REJECTIONS IN 5000 REPLICATIONS 
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CHI SQUARED TEST 
2 C PTS 3C PTS 50 PTS 





10* 


389 


439 


490 


CCNTFOL 


5* 


221 


251 


221 




1* 


38 


55 


45 


f- IX N CRM 


1C? 


1568 


2228 


3377 


VS 


5? 


1177 


1736 


2832 


N ( C,1 ) 


1? 


507 


551 


1756 







K0LMOG- SMR 


TEST 








20 PTS 


30 PTS 


50 PT 




10* 


535 


489 


490 


CCNTRCL 


5? 


258 


253 


245 




1* 


55 


54 


35 


VIXNCRM 


10* 


2467 


2194 


41C5 


VS 


5? 


1865 


2553 


3659 


N <C,1) 


1 % 


857 


1438 


2511 



FCUTZ FN TEST 

20 PTS 30 PTS 50 PTS 





10? 


502 


456 


495 


CCNTRGL 


5? 


2 52 


250 


235 




1? 


62 


51 




MXNCFM 


10? 


183 8 


2156 


2665 


VS 


5? 


1220 


1551 


1978 


N (0,1) 


1? 


486 


701 


975 



’CONTROL' IS THE TEST CF 



UNI FORM ( 0 »1 ) VS UN IFOR MO 1 1 ) 



•KI>NCRf» < IS TFE MIXED 
C. EC*N («5» 1) + 



NORMAL DISTRIBUTION 
0.2 0*N( C,4 ) 



NL'FBEF CF HYPOTHESIS REJECTIONS IN 5000 REPLICATIONS 







CHI SQUARED 


TEST 








20 PTS 


3 C D T S 


50 PTS 




10? 


431 


413 


5 2 6 


CCNTFCL 


5? 


240 


200 


245 




1? 


40 


42 


53 


PEARSON 


10? 


465 


478 


626 


VS 


5? 


258 


242 


302 


N (Ctl) 


1? 


46 


59 


71 







KOLMOG 


-SHIR TEST 








20 PTS 


3 C PTS 


5 0 PTS 




10? 


528 


478 


525 


CONTROL 


5? 


275 


242 


283 




1? 


45 


45 


54 


FEAR SCN 


10? 


623 


571 


673 


VS 


5? 


339 


310 


371 


N (C,l) 


1? 


66 


70 


55 






FOUTZ 


FN TEST 








20 PTS 


3C °TS 


50 FTS 




1C? 


512 


45 8 


5 3 C 


CONTROL 


5? 


248 


207 


263 




1? 


56 


48 


56 


PEARSON 


1C? 


481 


428 


457 


VS 


5? 


241 


196 


234 


N (0,1) 


1? 


60 


36 


57 



•CONTROL* IS IKE TEST OF UNIFORM<0, 1) VS UNIFORMS, 1) 

•PEARSON' IS THE PEARSON DISTRIBUTION 
BETA 1 = 0.0 , BETA 2 = 2.30 



MMBEP OF HYP C THE S I S REJECTIONS IN 5000 REPLICATIONS 
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CHI SOURED TEST 

20 PTS 3 C PTS 50 PTS 





10? 


435 


443 


541 


CCNTFCL 


5? 


237 


241 


247 




1? 


45 


52 


45 


PEAR SCN 


10? 


432 


445 


533 


VS 


5? 


243 


247 


241 


N (C,I) 


1? 


43 


55 


52 







KOLNCG-SNI R 


TEST 








20 PTS 


3 C PTS 


50 PTS 




1C? 


489 


528 


489 


CCNTFCL 


5? 


238 


252 


245 




1 % 


55 


50 


64 


FEAPSCN 


10 % 


456 


543 


503 


VS 


5 % 


242 


274 


262 


N (C,l) 


1 % 


61 


54 


61 



FOUTZ FN TEST 

20 PTS 3 C PTS 50 PTS 





10? 


536 


467 


487 


CCNTFCL 


5? 


275 


22 7 


240 




1? 


61 


34 


43 


FEAPSCN 


10? 


523 


441 


467 


VS 


5 % 


265 


217 


228 


N (C,l) 


1 % 


58 


33 


43 



•CCNTFCL’ IS The TEST OF UMpCRM(OtI) VS UNIFORM(Otl) 

'PEARSCN' IS THE PEARSON DISTRIBUTION 
BETA 1 = 0.0 » BETA 2 = 2.80 



M.NBEP CF HVPC THESIS REJECTIONS IN 50C0 REPLICATIONS 
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CHI SQUARED TEST 







20 FTS 


3 C p T S 


50 PTS 




1 c-s 


4C8 


390 


453 


CCNTFCL 


5? 


230 


214 


213 




1? 


35 


50 


52 


FEAPSCN 


less 


756 


992 


1576 


VS 


5? 


495 


634 


973 


N <C,1) 


1? 


121 


185 


343 







KOLNCG-SMI R 


TEST 








20 PTS 


3 C p TS 


50 FTS 




10? 


492 


515 


520 


CCNTFCL 


5? 


252 


255 


264 




1? 


55 


49 


52 


FEAPSCN 


10? 


942 


1161 


1498 


VS 


5? 


5 46 


693 


944 


N <C,1) 


1? 


158 


192 


294 







FOUTZ FN 


TEST 








20 PTS 


3 C PTS 


50 PTS 




10? 


5 05 


478 


5C4 


CCNTFCL 


5? 


2 23 


261 


2 e c 




1? 


50 


56 


48 


FEAPSCN 


10? 


8 53 


1400 


2487 


VS 


5? 


469 


884 


1785 


N (0,1) 


1? 


129 


311 


799 



•CCNTFCL * IS THE TEST CF UNIFCRM(0,1) VS U-NI FORM ( 0 , 1 ) 

•PEARSCN* IS THE PEARSON DISTRIBUTION 
BETA 1 = O.Ol, BETA 2 = 1.75 



NLME E F CF HYPOTHESIS REJECTIONS IN 5000 REPLICATIONS 
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CHI SQUARED TEST 

20 PTS 3C °T S 50 PTS 





10? 


428 


410 


452 


CCNTFGL 


5? 


239 


205 


2 a 




It 


40 


43 


52 


FEAFSCN 


10 % 


614 


744 


1082 


VS 


5% 


375 


43 5 


610 


N (C,l) 


1% 


93 


120 


153 







KOLMOG-SNIR 


TEST 








20 PTS 


3 C PTS 


50 PTS 




10? 


543 


470 


506 


CCNTFGL 


5? 


250 


229 


274 




1? 


41 


50 


48 


FEAFSCN 


10? 


837 


857 


1106 


VS 


5? 


482 


479 


64 8 


N <C,1) 


1% 


121 


131 


184 



FOUTZ FN TEST 



20 PTS 3C PTS 50 PTS 





1C? 


5 11 


474 


515 


CCNTFCL 


5? 


240 


214 


259 




1? 


55 


49 


54 


FEAR SON 


1C? 


558 


726 


1288 


VS 


5? 


305 


405 


774 


N (C,l) 


1? 


82 


1C1 


233 



•CCNTFCL' IS THE TEST QP UN I PQRM( C , 1 ) VS UNIFCRM<0,1) 

•PEARSCK' IS th£ PEAFSON DISTRIBUTION 
BETA 1 = 0.01, BETA 2 = 1.90 



M.H6EF OF HYPOTHESIS REJECTIONS IN 5000 REPLICATICNS 
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CHI SQUARED TEST 



20 PTS 3 C PTS 50 PTS 



CCMRCl 


10? 


403 


442 


510 


5? 


230 


237 


24! 




1? 


48 


41 


47 


FEAR SCN 


10% 


458 


556 


702 


vs 


5% 


267 


321 


366 


N (C,I) 


1% 


51 


74 


56 







KOLHOG-SVIR 


TEST 








20 PTS 


3 C PTS 


50 PTS 




10% 


497 


485 


5C2 


CCNTFCL 


5% 


240 


272 


258 




1% 


55 


50 


65 


PEAR SCN 


10% 


581 


650 


715 


VS 


5% 


3C9 


2 29 


365 


N (C,i) 


1% 


74 


73 


54 



FOUTZ FN TEST 



20 FT S 30 PTS 50 PTS 





10% 


525 


479 


546 


CCNTFCL 


5% 


241 


241 


285 




1% 


62 


56 


54 


PEARSCN 


10% 


534 


554 


666 


V S 


5% 


265 


283 


353 


N (C,l) 


1% 


68 


66 


81 



•CCNTFCL' IS THE TEST CF UNIF0RM(0,i) VS UNIFORMS, 1) 

•P EARSCN ' IS THE PEARSON DIS T RI8UTI0N 
BETA 1 = 0.25, BETA 2 = 3.20 



M.MEER OF HYPOTHESIS REJECTIONS IN 5000 REPLICATIONS 
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CHI SQUARED TEST 

20 PTS 3 C PTS 50 PTS 





10? 


373 


419 


502 


CONTROL 


5? 


22C 


230 


211 




1% 


55 


49 


24 


FEAPSCN 


10% 


739 


1C30 


1672 


VS 


5% 


476 


664 


1099 


N (C,l) 


\% 


118 


231 


421 







KOLHOG-SNIR 


TEST 








20 PTS 


3 C PTS 


50 PTS 




10% 


5 16 


483 


475 


CCNTFCL 


5% 


253 


227 


221 




1% 


48 


51 


38 


FEAPSCN 


10% 


790 


896 


1125 


VS 


5% 


474 


528 


735 


b <C,1 > 


1% 


121 


147 


224 



FOUTZ FN TEST 

20 PTS 3C PTS 50 PT$ 





1C? 


499 


481 


477 


CCNTFCL 


5% 


227 


238 


252 




1% 


46 


49 


5C 


FEAPSCN 


10% 


868 


1090 


1516 


VS 


5% 


471 


660 


965 


n icai 


1% 


121 


185 


344 



‘CCNTFCL' IS THE TEST OF UNI FQRM{ C ♦ 1) . VS UNIFORMS, I) 

' PEARSCN • IS THE PEARSON DISTRIBUTION 
EE T A 1 = 0 .5 0 1 SETA 2 = 3.00 



NLNSEF OF HYPOTHESIS REJECTIONS IN 50C0 REPLICATIONS 
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CHI SQUARED 


TEST 








20 PTS 


3 C PTS 


50 PTS 




1C? 


4 C 1 


428 


522 


CCNTFCL 


5? 


236 


225 


239 




1? 


43 


39 


54 


FEAPSCN 


10? 


2440 


3386 


46 54 


vs 


5? 


1960 


2820 


43 2 C 


N (C,l) 


1? 


931 


1632 


3267 







KOLMOG-SMIR 


TEST 








20 PTS 


3C PTS 


50 FT 




1C? 


482 


513 


509 


CCNTFCL 


5? 


239 


254 


245 




1? 


47 


52 


61 


FEAPSCN 


1C? 


11 55 


1459 


20C 1 


VS 


5? 


730 


920 


1428 


N < C,1 ) 


1? 


216 


322 


6C9 



FQUTZ F N TEST 

20 FT S 30 PTS 50 PTS 





10? 


531 


492 


55C 


CCNTFCL 


5? 


250 


250 


277 




1? 


57 


51 


61 


FEAPSCN 


1C? 


2135 


2837 


392C 


VS 


5? 


1460 


2156 


3309 


N < C ,1 ) 


1? 


583 


1049 


21C4 



'CCNTFCL' IS THE TEST CF UNIFORMIO,!) VS UNIFORMIO, 1) 

' P E ARSCN ' IS THE PEAFSON DISTRIBUTION 
BETA 1 = l.OC, BETA 2 = 3.40 



NLNE £ F CF HVPCTI-ESIS REJECTIONS IN 5000 REPLICATIONS 
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CHI SCUiREO TEST 

2C PTS 3 C P T S 50 PTS 





10? 


411 


427 


475 


CCNTFCL 


5? 


226 


238 


227 




1? 


34 


49 


34 


FEARSCN 


1C? 


1446 


2151 


4512 


\ S 


5? 


1070 


2546 


4035 


N (C,l) 


1? 


4C3 


1368 


2895 







KOL^OG-SMIR 


TEST 








20 PTS 


3 C PTS 


50 PTS 




10? 


495 


495 


5C4 


CCNTFCL 


5? 


249 


260 


265 




1? 


55 


64 


6 C 


FEARSCN 


1C? 


1032 


1314 


1765 


VS 


5? 


645 


832 


1227 


N (C,l) 


1? 


211 


265 


453 







FCUTZ FN 


TEST 








20 PTS 


3 C PTS 


50 PTS 




10? 


482 


481 


485 


CCNTFCL 


5? 


251 


251 


2C7 




1? 


51 


45 


43 


FEAFSCN 


10? 


1578 


2161 


3115 


vs 


5? 


1009 


1493 


2354 


N (C,l) 


1? 


348 


633 


1254 



•CCNTFCL' IS The TEST CP UNIF0RM(0,1> VS UNIFORMS, 1) 

• PEAR SCN ' IS THE PEARSON DISTRIBUTION 
BET* 1 = 1.00, BETA 2 = 3.60 



NLPBEF OF HYPCThESIS REJECTIONS IN 5000 REPLICATIONS 
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CHI SQUARED TEST 







20 FT S 


3 C PTS 


50 PTS 




1C? 


397 


427 


4 6 5 


CCNTFCL 


5? 


2 32 


233 


224 




1? 


41 


60 


56 


FEAPSCN 


1C? 


1086 


2009 


3226 


VS 


5? 


764 


1463 


2542 


N (C ,1 J 


1? 


255 


651 


1362 







KCLNCG-SMI R 


TEST 








20 PTS 


3 C PTS 


50 PTS 




10? 


524 


501 


459 


CCNTFCL 


5? 


267 


250 


252 




1? 


51 


56 


62 


FE/FSCN 


10? 


911 


1150 


1526 


VS 


5? 


5 46 


732 


1022 


N (Ctl) 


1? 


171 


229 


355 







F0UT2 FN 


TEST 








20 PTS 


3 C PTS 


50 PTS 




1C? 


521 


460 


525 


CCNTFCL 


5? 


’53 


248 


260 




1? 


50 


63 


46 


FEAP SCN 


10? 


1261 


1656 


2425 


VS 


5? 


764 


1070 


17C6 


N (C.l) 


1? 


243 


387 


725 



'CCNTFCL' IS TEE TEST CF UMFCRM(0,I) VS UNIFCRN<0»1> 

* PEARSCN * IS THE PEARSON DISTRIBUTION 
BETA I = l.OC, BETA 2 = 3. SO 



MReEF CF HYPOTHESIS REJECTIONS IN 50CO REPLICATIONS 
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CHI SQUARED TFST 

20 PTS 30 PTS 50 PTS 





10* 


371 


431 


464 


CONTROL 


5 % 


206 


237 


227 




1 % 


27 


48 


45 


FEAPSCN 


10* 


2231 


3C39 


4476 


VS 


5 % 


1745 


2282 


3538 


N i C » 1 ) 


1 % 


763 


1242 


2652 







KOLMCG-SMIR 


TEST 








20 PTS 


3 C D TS 


50 PTS 




1C* 


498 


482 


517 


CCNTFOL 


5* 


250 


261 


261 




1* 


42 


48 


54 


FEAPSCN 


10* 


11 53 


1481 


2082 


VS 


5* 


753 


1005 


1447 


n (cai 


1? 


221 


353 


621 



FOUTZ FN TEST 

20 PTS 30 DTS 50 PTS 





10* 


45C 


498 


486 


CONTROL 


5* 


258 


255 


248 




1* 


44 


45 


61 


FEAPSCN 


1C* 


1872 


2587 


3548 


VS 


5* 


1220 


1920 


2894 


N <C ,1 ) 


l* 


469 


886 


1644 



*CCNT PCI’ IS THE TEST CF UNIFORM(Otl) VS UNIFORMS, 1) 

•FEARSCN* IS THE PEARSON DISTRIBUTION 
BETA 1 = 2.0C, BETA 2 = 5.50 



NLMEEP OF HYPOTHESIS REJECTIONS IN 5000 REPLICATIONS 
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